







t*- ’ 









. T^=r:t;#Ts*^5: 








••Jt-t'^' r--. iSf * jty^tm -^‘ 




w ; 4 ig® 




A^_T*'? - . /ij: 





i-^:^ 







^1 






- s' 


















rnw:mm 


HYDES - A GENERALIZEIJ JI^YBRID 
COMPUTER PROGRAM FOR JTODYj 


»^mfM 

M 


TURBOJET OR TURBOFAN^lNGINl^y^ 


— - ^ _.._ . . 


Clevelmd, QMQ 


NATIONAL AERONA 



















. “ ‘ -3 *? 








r^TT^^-iEi- , . -; ■•' ■■■ . •;. 




e yfW^ 




:mmM 


-j^ 


-■yyrr^j-., 













1. Report No. 

NASA TM X-3014 


2. Government Accession No. 


Title and Subtitle , 

HYDES - A GENERALIZED HYBRID COMPUTER PROGRAM FOH| 
STUDYING TURBOJET OR TURBOFAN ENGINE DYNAMICS 


7. Author(s) 

JohnR. Szuch 


9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


3. Recipient's Catalog No. 


5. Report Date 

APRIL 197** 


6. Performing Organization Code 


8. Performing Organization Report No. 

E-7316 


10. Work Unit No. 

501-24 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

This report describes HYDES, a hybrid computer program capable of simulating one- spool 
turbojet, two-spool turbojet, or two-spool turbofan engine dynamics. HYDES is also capable 
of simulating two- or three- stream turbofans with or without mixing of the exhaust streams. 

The program is intended to reduce the time required for implementing dynamic engine simu- 
lations. HYDES was developed for running on the Lewis Research Center's Electronic 
Associates (EAI) 690 Hybrid Computing System and satisfies the 16 384-word core- size and 
hybrid- interface limits of that machine. The program could be modified for running on other 
computing systems. The use of HYDES to simulate a single-spool turbojet and a two- spool, 
two-stream turbofan engine is demonstrated. The form of the required input data is shown and 
samples of output listings (teletype) and transient plots (x-y plotter) are provided. HYDES is 
shown to be capable of performing both steady-state design and off-design analyses and transient 
analyses. 


17, Key Words {Suggested by Author(s) ) 

Hybrid computer; Simulation; Engine; Fan 


18. Distribution Statement 

Unclassified - unlimited 


Propulsion; Turbojet; Dynamics; Analog; 
Turbofan; Transient; Digital; Generalized; 
Control; Integration; Performance 


:at.28 

19. Security Classif, (of this report) 

Unclassified 

20. Security Classif. {of this page) 
Unclassified 

21. No. of Pages 

134 

22. Price* 

$ir.75 


' For sale by the National Technical Information Service, Springfield, Virginia 22151 






CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 2 

ENGINE CONFIGURATIONS 3 

MATHEMATICAL MODEL lO 

Intercomponent Volumes 10 

Fans and Compressors 11 

Turbines 14 

Rotor Dynamics 16 

Nozzles 16 

Combustor and Ducts 18 

HYBRID COMPUTER PROGRAM 19 

Hybrid Computing System 19 

Scaling 19 

Engine Definition 22 

Digital Program 23 

Data input program 23 

Main program 23 

Subroutines 24 

Analog Program 28 

Input Data Preparation 29 

Digital coefficients and DAC initial conditions 31 

Component performance data 32 

Analog coefficients 40 

EXAMPLES 40 

Single-Spool Turbojet (Configuration G) 40 

Turbofan (Configuration E) 50 

CONCLUDING REMARKS 59 

APPENDIXES 

A - UNSCALED-VARIABLE SYMBOLS 61 

B - SUMMARY OF UNSCALED EQUATIONS 68 

C - SCALED-VARIABLE SYMBOLS 83 

D - DEFINITIONS OF DIGITAL AND ANALOG COEFFICIENTS 88 

iii 



Page 


E - COMPUTER LISTINGS OF DATA INPUT PROGRAM, MAIN 

PROGRAM, AND SUBROUTINES 106 

F - SUMMARY OF SCALED EQUATIONS SOLVED ON ANALOG COMPUTER . . 120 

G - COMPUTER LISTING OF TURBINE DATA CONVERSION PROGRAM 124 

H - TURBOFAN FUEL CONTROL SIMULATION 127 

REFERENCES 131 


iv 


HYDES - A GENERALIZED HYBRID COMPUTER PROGRAM FOR 
STUDYING TURBOJET OR TURBOFAN ENGINE DYNAMICS 
by John R. Szuch 
Lewis Research Center 

SUMMARY 

The selection of an efficient airframe- engine combination depends on the ability to 
analyze a broad range of engine types and sizes operating at both design and off- design 
conditions. Computer programs having the necessary steady-state calculation capabili- 
ties have previously been developed. 

In addition to their use in steady-state studies, analyses of engine dynamics and 
control are also important in the selection of a suitable propulsion system. This report 
describes HYDES, a hybrid (analog- digital) computer program capable of simulating 
one-spool turbojet, two-spool turbojet, or two-spool turbofan engine dynamics. The 
program is also capable of simulating two- or three-stream turbofans with or without 
mixing of the exhaust streams. HYDES is intended to reduce the time required to im- 
plement dynamic engine simulations. 

HYDES was developed for running on the Lewis Research Center’s Electronic Asso- 
ciates (EAI) 690 Hybrid Computing System and satisfies the core- size (16 384 words) 
and hybrid-interface limits of that machine. The documentation of the program should 
allow the user to implement the program on another machine. The techniques employed 
and the resultant set of equations should be applicable in other generalized or specific 

engine simulations (analog, digital, or hybrid). 

For the existing program, procedures are described for generating the required in- 
put data, specifying the engine configuration of interest, and operating the program. 

The use of HYDES to simulate selected turbojet and turbofan engines is demonstrated. 
The form of the required input data for these engines together with samples of output 
listings (teletype) and transient plots (x-y plotter) for each example are provided. 
HYDES does not provide for engine control but does accept, as user-supplied inputs, 
fuel flow rates and nozzle areas. The simulation of simplified fuel control systems on 
the analog computer is discussed, and their use with the hybrid program is demon- 
strated. 



INTRODUCTION 


The selection of an efficient airframe- engine combination depends on the ability to 
analyze a broad range of engine types and sizes operating at both design and off-design 
conditions. Computer programs having the necessary steady-state calculation capabil- 
ities have previously been developed. The SMOTE code, discussed in references 1 and 
2, provides steady- state design and off- design calculation capability for both existing 
and theoretical turbofan engines. Theoretical engines are simulated by scaling compo- 
nent performance from existing engines to the design conditions of the theoretical 
engine. GENENG (ref. 3) extends the same techniques to handle turbojet engines as 
well as turbofans. GENENG II (ref. 4) was derived from GENENG and adds the capa- 
bility of studying two- or three- spool turbofan engines having as many as three nozzles 
(airstreams). 

In addition to their use in steady-state design- and off-design-point studies, anal- 
yses of engine dynamics and control are also important in the selection of a suitable 
propulsion system. For example, the use of turbofan engines as lift units for V/STOL 
aircraft (ref. 5) poses a number of engine control problems. At low flight speeds, the 
lift system must provide the fast thrust response needed for aircraft attitude control. 
The propilsion system must also be capable of correcting for upsetting moments caused 
by the loss of a lift engine. The required rapid engine accelerations must be accom- 
plished without exceeding turbine-temperature, rotor-speed, and compressor- stall 
limits. 

This report describes HYDES, a hybrid (analog-digital) computer program capable 
of simulating one- spool turbojet, two- spool turbojet, or two- spool turbofan engine 
dynamics. The program can easily be modified to handle the one-spool turbofan case. 
HYDES is capable of simulating two- or three-stream turbofans with or without mixing 
of the exhaust streams. The program is structured so as to allow the simulation of a 
wide range of engine sizes (as well as types) without changing the basic program. The 
hybrid computer was used because it combines the precision and logic capabilities of 
the digital computer with the integration and output capabilities of the analog computer. 
Steady-state design-point data (pressures, temperatures, etc. ), generated by a program 
such as GENENG II, and the associated component performance maps are required to 
generate the necessary input data for the hybrid program. After the simulation has 
"settled-out" at the design point, a change in one of the input variables (such as fuel 
flow) will result in a transient excursion to a new steady-state operating point. Thus, 
the program is capable of steady-state design- and off-design-point operation, as well 
as x»ssessing the desired transient capability. 

HYDES does not provide for the simulation of control system dynamics. The dig- 
ital portion of the program accepts values for fuel flow rates and nozzle area (or areas) 
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from the analog portion of the program. In general, the function generation and arith- 
metic operations required in simulating an engine control system could be performed 
on either an analog or digital computer. But hybrid interface limitations dictated that 
control system d 5 mamics (when required) be simulated on the analog computer. To de- 
monstrate the transient operation of the program for a selected engine, a fuel control 
system was assumed. The implementation of the fuel control and its relationship to 
the basic HYDES program are illustrated in this report. 

HYDES was developed for running on the Lewis Research Center's Electronic Asso- 
ciates (EAI) 690 Hybrid Computing System. The structure of the program is very much 
influenced by the digital core size (16 384 words) and the hybrid-interface capability 
(24 analog-to digital converters and 24 digital-to-analog converters) of that machine. 

The documentation of the program, together with the associated digital computer soft- 
ware (available from the author upon request), should allow the direct use of the program 
on another EAI 690 computer. The simulation techniques and resultant equations, pre- 
sented in this report, should also serve as a guide in the development of both generalized 
and specific turbojet and turbofan engine simulations for use with different computing 
systems. 


ENGINE CONFIGURATIONS 

HYDES can be used to simulate a number of different engine configurations. These 
configurations are referred to as configurations A to H. A discussion of each follows. 

A schematic representation of configuration A is shown in figure 1. Configuration 
A represents an unmixed, two- spool, three- stream turbofan engine with separate per- 
formance maps for the fan- hub (low-pressure compressor) and fan-tip sections. All 
other configurations can be considered as variations of this configuration. The fan is 
driven by a separate low-pressure turbine. The fan-hub flow Wfj discharges into a 
high-pr;^ssure- compressor-inlet volume V 2 which supplies (1) bleed flow for control 
purposes (2) flow to the second stream Wy2> (3) bleed flow to the third stream 

Wbig, and (4) flow to the core compressor w^. All unscaled-variable symbols are de- 
fined in appendix A. The core (high-pressure) compressor, which is driven by the 
high-pressure turbine, discharges into a combustor-inlet volume Vg which supplies 
(1) overboard bleed flow (2) cooling flow for the high-pressure turbine Wj^^j, 

(3) cooling flow for the low-pressure turbine Wj^j 2 ’ airflow to the combustor w^. 

The combustor airflow reacts in high-pressure-turbine-inlet volume with the 
injected fuel flow Wp. The high-pressure-turbine flow W|.j^ discharges into low- 
pressure-turbine-inlet volume Vg, where it is diluted by the high-pressure-turbine 
cooling flow Similarly, the low-pressure-turbine flow w ^2 discharges into 
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Figure 1. - Schematic representation of configuration A. 

core-duct volume Vg where it is diluted by the low-pressure-turbine cooling flow 
Wbi2‘ All turbine cooling bleeds are assumed to enter downstream of the rotor (or 
rotors), thus contributing no work. Volume Vg supplies flow to the core nozzle 
The second- stream flow w^2 discharges into the second- stream-duct volume V2 2 
which, in turn, supplies flow to the second- stream nozzle w^^^. The fan-tip flow w^2 
discharges into the bypass-fan-discharge volume Vg 2» where it is mixed with the 
inter stream bleed flow w^jg- Volume Vg g supplies flow to the third stream w^g. 

This flow discharges into the bypass- fan- duct volume Vg g which, in turn, supplies 
flow to the third- stream nozzle Constants and Kjg may be set to nonunity 

values to represent a gear-driven low-pressure compressor and fan, respectively. If 
the low-pressure compressor represents the fan hub, and K^g are equal. 

A schematic representation of configuration B is shown in figure 2 . This configu- 
ration is identical to configuration A except for the mixing of the core and third- stream 
flows in the mixing volume This volume then supplies flow to the mixed-flow nozzle 
Wn3‘ The second- stream nozzle flow w^^ discharges to the atmosphere at station 8 . 

A schematic representation of configuration C is shown in figure 3 . Configuration C 
represents a two- spool, two- stream turbofan engine with the entire fan represented by 
one set of performance maps. Thrust is generated by the separate exhaust flows w^^j 
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Figure 3. - Schematic representation of configuration C. 


















and This configuration is obtained by eliminating the fan-tip and third-stream 

calculations from configuration A (fig. 1). The mixed version of the two-spool, two- 
stream turbofan engine is shown, schematically, in figure 4 and is referred to as config- 
uration D. 

A schematic representation of configuration E is shown in figure 5. Configuration E 
represents a two- spool, two- stream turbofan engine with separate performance maps 
for the fan hub and tip sections. This configuration is formed by eliminating the second- 
stream calculations from configuration A. The mixed version of configuration E is 

Ps 



Figure 4. - Schematic representation of configuration D. 



6 




referred to as configuration F and is shown, schematically, in figure 6. 

A conventional one- spool turbojet engine is referred to as configuration G and is 
shown, schematically, in figure 7. The one-spool turbojet configuration is formed from 
configuration C (fig. 3) (1) by eliminating the fan and second-stream calculations, (2) by 
equating the conditions in volume V 2 j with the inlet conditions at station 2, (3) by 
eliminating the low-pressure- turbine calculations, and (4) by equating the conditions in 
volumes Vg and Vg. 

A two- spool turbojet is referred to as configuration H and is shown, schematically, 
in figure 8. This configuration is formed from configuration E (fig. 5) by eliminating 
the fan- tip and third- stream calculations. 

For configurations E and F, an option (-2) is provided in the HYDES program to 
allow the simulation of the fan tip without simulation of the fan hub. An example of this 
type of engine would be an aft-fan engine with a single core (high-pressure) compressor. 
For this option, the program does provide for supercharging of the high-pressure com- 
pressor (if required). A supercharger pressure ratio is computed as a linear function 
of the fan corrected speed. The supercharger torque is computed by using the high- 
pressure compressor flow (volume V 2 j is assumed to be negligible). If super- 
charging is not required, the appropriate coefficients are set to give a supercharger 
pressure ratio of 1.0. When a gear-driven aft- fan is being simulated, the constant 
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Kj 2 is set equal to the gear ratio. 

Another option (-1) is provided for those configurations requiring separate perfor- 
mance maps for the fan hub and tip sections (A, B, E, and F). For these configurations, 
the fan-tip performance data may be in terms of either the total fan flow (w^j + w^g) or 
the fan- tip flow w^g- Table I summarizes the configurations and options available with 
the HYDES program. The simultaneous use of the -1 and -2 options is referred to as 
the -3 option. 


TABLE I. - ENGINE CONFIGURATIONS AND OPTIONS FOR HYDES PROGRAM 


Configuration 

Description 

A 

Two-spool, three- stream, unmixed turbofan; separate maps for fan 
hub and tip sections 

A-1 

Same as configuration A with fan-tip performance in terms of total 
fan flow 

B 

Two-spool, three- stream, mixed turbofan; mixing of first and third 
streams; separate maps for fan hub and tip sections 

B-1 

Same as configuration B with fan-tip performance in terms of total 
fan flow 

C 

Two-spool, two-stream, unmixed turbofan; entire fan described by 
same map 

D 

Two-spool, two-stream, mixed turbofan; entire fan described by the 
same map 

E 

Two- spool, two- stream, unmixed turbofan; separate maps for fan hub 
and tip sections 

E-1 

Same as configuration E with fan-tip performance in terms of total 
fan flow 

E-2 

Same as configuration E with fan hub performance simplified for super- 
charger or aft- fan 

E-3 

Same as configuration E with fan-tip performance simplified for super- 
charger or aft-fan 

F 

Two- spool, two- stream, mixed turbofan; separate maps for fan hub 
and tip sections 

F-1 

Same as configuration F with fan-tip performance in terms of total 
fan flow 

F-2 

Same as configuration F with fan-hub performance simplified for super- 
charger or aft-fan 

F-3 

Same as configuration F with fan-tip performance simplified for 
supercharger or aft-fan 

G 

One-spool turbojet 

H 

Two- spool turbojet 
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MATHEMATICAL MODEL 


The first step in developing any simulation is the formulation of a mathematical 
model. This model, in equation form, represents the functional relations that exist 
between system variables. In the case of turbojet and turbofan engine systems, these 
variables are pressures, temperatures, flow rates, rotor speeds, and so forth. 


Inter component Volumes 

Pressures and temperatures are computed in each of the intercomponent volumes 
shown in figures 1 to 8. In these volumes, storage of mass and energy occurs (ref. 6). 
Modified forms of the continuity and energy equations written for each volume are 



( 2 ) 


A summary of all equations written for specific components is given in appendix B. 

With the results from equations (1) and (2), the pressure in the volume can be computed 
from the ideal-gas law: 


P=:-WT (3) 

V 

The thermodynamic properties of air and fuel-air mixtures (c^, c^, 7 , h) are 
calculated by considering variable specific heats and no dissociation. The air and 
fu 0 l_air property tables of reference 7 were curve-fit by the authors of references 1 
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and 2. Those curve-fits were further simplified and used in the HYDES program. For 
each intercomponent volume in a particular engine configuration, the following gas prop- 
erties are calculated: 


II 

(4) 

R = fg (f/a) 

(5) 

R 


c = c - — 
V P J 

(6) 

c 


II 

(7) 

"-V 


= fg (T,f/a) 

(8) 


where f/a is the local fuel-air ratio. 

While the gas constant R is, in general, a variable when mixtures of gases are 
considered, it was determined that the sensitivity of R to the fuel-air ratios expected 
in this type of simulation could be neglected. Therefore, the gas constant of air R^ 
is used in equation (3). The use of a constant value of R also prevents the occurrence 
of algebraic loops, which require iterative solutions. 


Fans and Compressors 

Fans and compressors can be modeled by a number of known techniques. One meth- 
od is to represent multistage compressors with individual stage models (i. e. , compute 
pressure and temperature rises across each stage). This technique is referred to as 
stage- stacking (ref. 8), but it requires a large computing facility when used in a total 
engine simulation. For the HYDES program, fans and compressors are represented 
by overall performance maps. This technique does not consider interstage gas dynam- 
ics. The effects of interstage bleeds or variable geometry, if they exist, must be re- 
flected in the overall performance data. Figure 9 shows the form of the fan and com- 
pressor maps used in this program. Figure 9(a) shows a corrected flow parameter fpfc 
plotted as a function of two variables - a pressure ratio pr and a corrected speed pa- 
rameter fcnp. When performance data are not available for a particular fan or com- 
pressor (as in the case of theoretical engines), it is necessary to scale available data 
to the design point of the theoretical component. Scale factor (WACF, ETACF, PRCF) 
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Map pressure ratio, pr 
(b) Efficiency map. 

Figure 9. - Form of fan and compressor perfor- 
mance maps used in hybrid simulation. 



are used to accomplish the performance scaling within the HYDES program. The valid- 
ity of the scaled data decreases as the scale factors deviate from unity. Figure 9(b) 
shows fan or compressor efficiency eff plotted as a function of the same two independ- 
ent variables. The following equations are used to compute fan or compressor flow, 
discharge enthalpy, and torque from specified inlet conditions, rotor speed, and back- 
pressure; 


pr 





in 


PRCF 


+ 1 


fcnp = 



( 9 ) 


( 10 ) 


12 


fpfc = f^(pr, fcnp) 


( 11 ) 


w = 


(WACF)(fp£c) 


P. 

in 




eff = fgCpr, fcnp) 


T] = (ETACF)(eff) 
7-1 




7TN 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 


The fan or compressor discharge temperature T^^^ and the corresponding en- 
thalpy h^^^ represent the inlet conditions to the downstream volume (eq. (2)). If the 
fan or compressor is the only component feeding the volume, the discharge conditions 
^out ^out equal (in steady state) the temperature and enthalpy, respectively, 
in the doAvnstream volume. 

The solution of equation (15) requires a knowledge of the "average" thermodynamic 
properties in the fan or compressor. Since these properties are functions of temper- 
ature, they vary throughout the component. For this reason, a temperature interpo- 
lation constant is adjusted for each fan or compressor to match available steady - 
state cycle data for the engine being simulated. That is. 


T=^T.^.(1-^T^^^ 


(18) 
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= f^(T, f/a) 

(19) 

R = f2(f/a) 

(20) 

c =c -H 
^ P J 

(21) 

c 

y 

(22) 

c 

V 



The calculated value for y should agree with the specific-heat ratio used to define the 
fan or compressor efficiency rj. The temperature in the discharge volume T^^^. is 
used in calculating T to avoid the occurrence of algebraic loops associated with the 
use of 


Turbines 

As in the case of fans and compressors, the most direct approach to modeling 
multistage turbines would be to apply stage -stacking techniques. However, individual 
stage performance data are usually not available. It is therefore necessary to repre- 
sent turbines by overall performance maps. Figure 10 shows the form of the turbine 
maps used in the HYDES program. Figure 10(a) shows a turbine flow parameter fpt 
plotted as a function of pressure ratio pr and a turbine speed parameter tnp. 

Figure 10(b) shows a turbine enthalpy drop (work) parameter hpt plotted as a function 
of the same two independent variables. As in the case of fans and compressors, scale 
factors (CNCF, DHCF, TFCF) are used to scale available performance data to the design 
point of theoretical turbines. The following equations are used to compute the turbine 
flow, discharge enthalpy, and torque from specified inlet conditions, rotor speed, and 
back- pressure: 


pr 



(23) 
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(b) Enthalpy drop (work) map. 

Figure 10. - Form of turbine performance maps 
used in hybrid simulation. 


tnp = <CNCF) (N) 
1 in 

(24) 

= fg(pr, tnp) 

(25) 

hpt = f,^(pr, tnp) 

(26) 

(CNCF)(^t)(P^ )(N) 
(TFCF)(T,„) 

(27) 


Ah 


(DHCF)(hpt)(N)Y/T~ (CNCF) 
1000 


(28) 


out 


= h. - 
in 


Ah 


(29) 
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_ 30J(4h)(w) 
7TN 


( 30 ) 


It should be noted that a factor of 1000 has been included in the definition of the enthalpy 
parameter hpt. 

In general, the discharge enthalpy for the turbine h^^^ does not equal (in steady - 
state) the enthalpy in the downstream volume h^^^ because of the mixing of the turbine 
flow with the turbine cooling bleed. 


Rotor Dynamics 

After the fan, compressor, and turbine torques are computed, the rotor speed 
(or speeds) is computed by using the conservation of angular momentum. That is, 

N=— f ALdt + N. (31) 

7TI yp ^ 

where AL denotes the difference between the driving turbine torque and the load torque 
(fans or compressors) and I represents the polar moment of inertia for the spool. 


Nozzles 

For all the engine configurations being considered (figs. 1 to 8), flows must be 
computed for each nozzle. All nozzles are assumed to be the convergent type, and the 
flow processes are assumed to be isentropic. For a specified inlet pressure, there 
exists, for each nozzle, a critical back -pressure (ref. 6) given by 

y 

P = P- (32 

\y + i/ 

If the back -pressure is higher than the critical pressure, the flow is subsonic at 
the nozzle throat and may be expressed as 
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1 


w = 





y-1 



( 33 ) 


If the back-pressure is lower than the critical pressure, the flow is sonic or "choked” 
at the nozzle throat. For this case, the flow is given by 


w 




( 34 ) 


For those nozzles where thrust calculations are required, the following equations are 
used; 



for subsonic flow, and 



for sonic flow. 

To minimize the digital computation time, the effects of varying Y and R in the 
nozzle flow and thrust calculations have been neglected in the HYDES program. A 
constant value of y, 1. 35, and the gas constant of air are used for all nozzle calcu- 
lations. 

The turbine cooling, control, and overboard bleed flows are treated in the same 
manner as the nozzle flows. However, the high pressure ratios between the compressor 
discharge volume Vg and the turbine discharge volumes and Vg and the ambient 

pressure allow the turbine cooling and overboard bleed flows to be computed by using 
equation (34) without testing the pressure ratio. 
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Combustor and Ducts 


For all engine configurations, total pressure losses are assumed in the combustor 
and in each of the bypass ducts (for turbofans )l Mach numbers in the combustor and 
ducts are assumed to be low enough to allow the assumption of incompressible flow. 
Flows are computed from 



where 01 is the flow resistance as determined from steady-state design-point data. 

The inlet and discharge pressures for the combustor and bypass ducts are computed 
by using equations (1) to (3). For the bypass ducts, no heat addition to the ducts is 
assumed. Therefore, the input enthalpy for the duct-discharge temperature calculation 
IS the enthalpy in the duct inlet volume. For the combustor, heat addition due to the 
burning of the injected fuel is assumed to take place in the combustor discharge volume 
V 4 . In the calculation of the combustor discharge temperature T^, the input energy 
to the volume is given by 


NN 

E 


w. . h. 
in, j in, ] 


= Wj,hb+n 


bWpHVF 


(38) 


where hj^ denotes the enthalpy of the combustor air. For this program, a constant 
combustor efficiency 77 ^ is assumed. An interpolation constant & is adjusted to 
satisfy steady-state cycle data. That is, 

-^b>T4 

(39) 

\ = h Pb. ((^>4] 

(40) 

li 

(41) 
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HYBRID COMPUTER PROGRAM 


Hybrid Computing System 

The HYDES program was developed for running on the Lewis Research Center's 
Electronic Associates (EAI) Model 690 Hybrid Computing System. The 690 System 
(ref. 9) consists of an EAI 640 Digital Computer, an EAI 693 Hybrid Interface Unit, 
and an EAI 680 Analog Computer. 

The basic digital computer has 16 384 words of core storage. Floating-point (real) 
numbers are represented by two 16-bit words. Scaled-fractions (numbers less than 1. 0) 
are represented by a single 16-bit word. Arithmetic operations can be performed by 
the digital computer for both floating-point and scaled- fraction numbers. For digital 
input and output, the hybrid computer uses a teletypewriter, a high-speed paper-tape 
reader, and a high-speed paper-tape punch. 

The interface unit provides the necessary communication between the analog and 
digital computers. Twenty-four analog-to-digital converters (ADC's) are used to 
transmit analog signals to the digital. Twenty-four digital-to-analog converters (DAC's) 
are used to transmit digital signals to the analog. The interface system also contains 
16 control lines and eight sense lines. The logical states of the control lines are set 
by the digital program and may be sensed on the analog. Similarly, the logical states 
of the sense lines are set by the analog and sensed by the digital. Eight sense switches 
can be positioned at the digital control console and tested by the digital program. 

The analog computer performs all the operations characteristic of analog machines 
(i. e. , summing, integration with respect to time, limiting, attenuation, multiplication, 
function generation, etc. ). The analog computer contains a total of 156 amplifiers, 

124 potentiometers, and 24 quarter-square multipliers. In addition, the analog contains 
16 comparators and 16 function relays which can be positioned by either the digital or 
analog computers. The use of peripheral equipment such as x-y plotters and strip- chart 
recorders allows continuous monitoring by the user of computed variables. 


Scaling 

To reduce the core requirements and computation time of the digital portion of the 
HYDES program, it was decided to use scaled fractions throughout the digital program. 
Therefore, all digital variables are scaled so as not to exceed unity during the program 
execution. For each variable x, a scale factor SF^^ is chosen so as to limit the scaled 
variable X = x/SF^ to the range -1<X< +1. Similarly, all analog variables are scaled 
so that no analog signal exceeds 1.0 computer unit (10 volts on the EAI 680 computer). 



While the choice of most scale factors is left to the user, certain variables have 
been prescaled to minimize core storage requirements. For example, all fuel-air 
ratios are scaled for a maximum of 0. 05. In addition, certain variables are assumed to 
have the same scale factors as other variables. Table n contains a list of prescaled 
variables, their respective scale factors, and the implied relations between scale 
factors. 

In addition to the previously described amplitude scaling, it is often necessary to 
time-scale the analog portion of a dynamic simulation. To allow the treatment of dig- 
ital outputs as continuous input signals to the analog, a sufficient number of cycles 
through the digital loop must occur for each cycle of the analog frequencies. In some 
cases, computational stability can only be achieved by decreasing the analog frequencies. 
A time-scale factor SF^., is selected such that the computer time t' equals SF^, t. For 

TABLE n. - PRESCALED VARIABLES FOR 


HYDES PROGRAM 


Variable, x 

Scale factor^, SF^ 

2’ ^2.1’ ^2.r ^3.2’ 


'^3.2’ '^2.2’ '^3. 5' ’^tl’ ' 

555.55 K (1000° R) 

"^£2 


^3- ^3 T, 

1111.1 K (2000° R) 

^4' Tb 

2777.8 K (5000° R) 

TS’ Ts 

2222.2 K (4000° R) 

T^, Tg, 

1388.9 K (2500^ R) 


SF^jSF^p^ 


2092. 2 J/kg-K (0. 5 Btu/lbm-^R) 

r 

2.0 


1.0 

^2 

1.0 

tr 

1.25 

f/a 

0.05 


^The following relations between scale factors are assumed: 
^^TRQT1"®^TRQC 

®^WDF1"^^WDC"®^WDB=®^WDT1"®^WDT2^®^WDN1 

®^WDF2^®^WDV3^^^ WDN2 

^^WDV2=®^WDN4 

SFpg=SFp4 

®^P32’®^P35 
SFps=SFpiN for nozzles 
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this definition of SF^.,, a value exceeding 1. 0 corresponds to a "slowing down" of the 
analog problem. 

The selection of a suitable time-scale factor requires the user to estimate the real- 
time frequencies associated with the engine dynamics. Analog frequencies can be esti- 
mated for each inter component volume by assuming a "resistive" termination down- 
stream. For this simplified first-order model, the time constant associated with each 
volume is given by 


r - W _ PV _ 1 

w RTw 0> 
a 


(42) 


where W is the stored mass in the volume, w is the total flow through the volume, 
and is the real-time analog frequency. The ratio of digital frequency to the scaled 
cutoff frequency associated with the volume is given by 


0) 


d 



U) 

a 


(43) 


where t^ is the cycle time for the digital program. The cycle time for the existing 
digital program varies depending on the engine conf^ration but is between 23 and 
44 milliseconds. Experience has shown that a minimum frequency ratio of 50:1 at the 
design point results in computational stability. A smaller ratio might prove satisfactory 
for a particular engine simulation, however. The achievement of this minimum ratio 
usually requires a combination of time scaling (SF^.,) and increasing of the smaller inter- 
component volumes V . When increasing the volume above the actual value, care should 
be taken to keep the decreased cutoff frequency above the range of frequencies of 
interest. 

For solution on the analog or digital computer, the equations given in appendix B 
must be rewritten in terms of scaled variables. Appendix C contains the definitions 
of the scaled (computer) variables. For the scaled equations solved on the digital com- 
puter, the scale factors and engine parameters are combined to form digital coeffi- 
cients SC(i). Similarly, the scale factors and engine parameters appearing in the 
analog-solved equations are combined to form analog coefficients C(i). Appendix D con- 
tains the definitions of the digital and analog coefficients. 
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Engine Definition 


The use of the HYDES program to simulate any of the configurations described in 
table I requires the user to supply to the computer information which defines the engine 
to be studied. The engine configuration is specified by a set of logical variables. These 
variables are either "TRUE" or "FALSE." The means of setting these variables is dis- 
cussed in the section Input Data Preparation. Table HI lists the combinations of these 
variables which define the 16 configurations (with options) previously discussed. Combi- 
nations of the variables which do not appear in table m may lead to fallacious results 
and should be avoided. 

The logical variable HBPR="TRUE" denotes a fan requiring separate performance 
maps for the fan hub and tip sections (configurations A, B, E, and F). The variable 
MIX="TRUE” denotes mixed turbofans (configurations B, D, and F). The variable 
STRM3="TRUE" denotes three-stream turbofans (configurations A and B). 
TURBJ1="TRUE" and TURBJ2="TRUE" denote one-spool and two-spool turbojets, 
respectively (configurations G and H). For the HBPR="TRUE" configurations, 
SPLIT="TRUE" denotes the fan-tip performance data expressed in terms of total fan 
flow (-1 option). For the case when HBPR="TRUE" and STRM3="FALSE", SUPER= 
"TRUE" denotes the simplified fan-hub performance map previously discussed (-2 


TABLE m. - LOGICAL VARIABLES DENOTING ENGINE CONFIGURATIONS 


AND OPTIONS FOR HYDES PROGRAM® 


Configuration 

Logical variable 

HBPR 

MIX 

STRM3 

TURBJl 

TURBJ2 

SUPER 

SPLIT 

A 

1 


F 


T 

F 


F 

I 


F 

A-1 



F 










T 

B 



1 










F 

B‘l 



1 










T 

C 

I 


F 


F 







F 

D 

F 

1 










F 

E 

T 

F 










F 

E-1 













T 

E-2 











T 

F 

E-3 











T 

T 

F 



T 







F 

F 

F-1 











F 

T 

F-2 











T 

F 

F-3 












r 

T 

G 

F 

F 



T 



F 

F 

H 

] 


F 

1 


F 

1 


F 

F 


®T denotes "TRUE” and F denotes "FALSE.” 
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option). 

To simulate a particular engine, the user must supply the values for the analog and 
digital coefficients defined in appendix D. The means of inputting these coefficients to 
the computer is discussed in the section Input Data Preparation. The digital coefficients 
are assumed to be scale fractions (less than unity). If a calculated value exceeds unity 
for the selected scale factors and/or engine parameters, the coefficient must be rede- 
fined and the corresponding equation (or equations) modified accordingly. Similarly, 
the analog coefficients defined in appendix D represent potentiometer settings which 
should be less than 1. 0. Redefinition of an analog coefficient requires a corresponding 
change in an amplifier gain on the analog computer. The definitions of the analog and 
digital coefficients listed in appendix D have been chosen so as to minimize the changes 
required to simulate a wide range of engine types and sizes. 


Digital Program 

The digital portion of the HYDES program is used, primarily, to compute the time 
derivatives associated with the engine system's dynamics. The bracketed terms in 
equations (1) and (2) and the torque differentials associated with equation (31) are com- 
puted on the digital and transmitted to the analog by means of the DAC's. The analog 
computer is then used to compute the pressures, temperatures, and rotor speeds which 
are required as input to the digital program and are input through the ADC's. 

Because of the large size of the simulation, the digital program was divided into 
two parts - a data input program and the main program. 

D ata input program . - The data input program reads and stores (1) digital coeffi- 
cients SC(i), (2) initial conditions for the DAC's AI(i), (3) scale factors for the compon- 
ent map data, and (4) unsealed component performance data. The program proceeds to 
scale the component map data and stores the SC(i)'s, AI(i)'s, and scaled map data in 
COMMON blocks shared by the main program and associated subprograms. The pro- 
gram then tests the position of user-set sense switches to determine the states of the 
logical variables HBPR, MIX, etc. (see table III). Based on the specified configuration, 
the scaled component map data are stored in the proper storage locations. 

The data input program (and associated library routines) requires 7367 words of 
core storage. Included in this total are 2214 words of storage used for the COMMON 
data. Appendix E contains a FORTRAN listing of the data input program as written for 
the Lewis Research Center's hybrid system. (The reading-in of the scaled-fraction 

SC(i) and AI(i) data requires the use of a special S format (ref. 10) in statement 5 of 
the listing. ) 

Mam program . - The main digital program is used to perform all the algebraic 
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computation and function generation required to compute the temperature, stored-mass, 
and rotor-speed derivatives. The main program also tests sense switches to determine 

the engine configuration and proper branching within the program. 

To reduce the core storage requirements, the main program utilizes subroutines to 
perform repeated operations. For example, a subroutine FNCMP is used to solve the 
scaled versions of equations (9) to (17) for each of the fans and/or compressors. The 
following section describes each of the subroutines utilized by the mam program. The 
main digital program (and associated subroutines and library routines) requires 1 
words of core storage. Included in this total are the 2214 words of storage used for 


COMMON data. 

The main program accepts current values of pressures, temperatures, rotor 
speeds nozzle areas, and fuel flow as input from the analog computer. Updated values 
for temperature derivatives, stored- mass derivatives, torque differentials, compressor 
map variables, and thrust are output to the analog through the DAC’s. The user, by 
depressing a sense switch, can direct the program to output current values of selected 
engine variables on the teletype. The program also tests for scaled-fraction overflows 
and wUl output an appropriate message at the teletype if an overflow occurs. Appendix 
E contains a FORTRAN listing of the main program as written for the Lewis hybrid 
system (The use of scaled fractions in the main program requires the suffix S on 
fractional constants and a special routine SSQRT to perform the square root operation 


(ref. 10)). 

Subroutines - The various subroutines called by the main digital program are 
listed here anTthe function of each Is described. FORTRAN listings of the subroutines 

are included in appendix E. 

PROCOM- Subroutine PROCOM calculates thermodynamic gas properties either for 
air or for a JP4-air mixture using curve-fits of the tables of reference 7. Inputs to the 
subroutine are the scaled-fraction temperature (rescaled for a maximum of 277. 8 K 
(5000° R) prior to being input) and the scaled-fraction fuel-air ratio (scaled for a maxi- 
mum of 0. 05). The outputs of the subroutine are the scaled- fraction specific heats 
(scaled for a maximum of 0. 5), the specific-heat ratio (scaled for a maximum of 2. 0), 
and enthalpy (scaled for a maximum of 3.487x10® J/kg (1500 Btu/lbm). A FORTRAN 

listing of subroutine PROCOM is given in appendix E. 

THAT; Subroutine THAT calculates the isentropic temperature rise parameter 
(AT/T. ).. for a fan or compressor. The subroutine is called by the main program 
when SUPER=”TRUE” and by FNCMP. The inputs to TRAT are the scaled- fraction 
pressure ratio (rescaled for a maximum of 15. 0 prior to being input) and the average 
specific-heat ratio (scaled for a maximum of 2.0). The output of the subroutine is the 
scaled-fraction temperature rise parameter (scaled for a maximum of 1.25). Subroutine 
TRAT makes use of calls to subroutines NMXTR and CLCTR to perform the necessary 
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calculations. A FORTRAN listing of subroutine TRAT is given in appendix E. 

NOZZL: Subroutine NOZZL computes those terms, sensitive to pressure ratio and 
specific-heat ratio, required in the calculation of flow and thrust for a convergent nozzle. 
The inputs to the subroutine are the scaled- fraction inlet and exit pressures and the 
ratio of the exit and inlet pressure scale factors. The outputs of the subroutine are the 
functions 



(both scaled for a maximum of 1. 0) and the static pressure at the nozzle throat P 

s 

(scaled the same as The specific -heat ratio r is assumed to be 1. 35 for these 

calculations. Subroutine NOZZL makes use of calls to subroutines NRMPN, CLCFN, 
and CLCWN to perform the necessary calculations. A FORTRAN listing of subroutine 
NOZZL is given in appendix E. 

FNCMP: Subroutine FNCMP calculates the flow, discharge enthalpy, and torque 
for a fan or compressor. The inputs to the subroutine are a component index N (N=l 
for the fein hub or low-pressure compressor, N=2 for the high-pressure compressor, 

N=3 for the fan tip) and the scaled- fraction inlet pressure, inlet temperature, inlet 
enthalpy, exit pressure, exit volume temperature, and fraction of design rotor speed 
(the latter scaled for a maximum of 2. 0). The outputs of the subroutine are the scaled- 
fraction flow, discharge enthalpy, torque, map pressure ratio, and map flow parameter. 
Subroutine FNCMP makes use of calls to subroutine MAPFUN to perform the radial 
interpolation of the map data stored in COMMON arrays. A FORTRAN listing of sub- 
routine FNCMP is given in appendix E. 

TURB; Subroutine TURB calculates the flow, discharge enthalpy, and torque for 
turbines. Subroutine TURB is similar to subroutine FNCMP in operation. The inputs 
to TURB are a component index N (N=l for the high-pressure turbine, N=2 for the low- 
pressure turbine) and the scaled- fraction inlet pressure, inlet temperature, inlet 
enthalpy, exit pressure, and fraction of design rotor speed (scaled for a maximum of 
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2.0). The outputs of the subroutine are the scaled- fraction flow, discharge enthalpy, 
and torque. The FORTRAN listing of subroutine TURB is given in appendix E. 

NRMPN: Subroutine NRMPN converts a scaled-fraction pressure ratio to a format 
suitable for use with the library-supplied function generation routines. Subroutine 
NRMPN is called by subroutine NOZZL and contains the pressure-ratio breakpoints 
used in the computation of the nozzle flow and thrust functions. The input to the sub- 
routine is the scaled-fraction pressure ratio (scaled for a maximum of 1. 0). The out- 
puts of the subroutine are the normalized pressure ratio and an integer ERRFLG used 
to signal an out-of-range input. A FORTRAN listing of NRMPN is given in appendix E. 

CLCFN: Subroutine CLCFN generates the function 



required for convergent nozzle thrust calculations. Subroutine CLCFN is called by sub- 
routine NOZZL and contains the function values (scaled for a maximum of 1. 0) corres- 
ponding to the breakpoints contained in NRMPN with V equal to 1. 35. The input to the 
subroutine is the normalized pressure ratio. The output of the subroutine is the scaled- 
fraction thrust function. A FORTRAN listing of CLCFN is given in appendix E. 

CLCWN: Subroutine CLCWN generates the function 



required for convergent nozzle flow calculations. CLCWN is called by subroutine NOZZL 
and contains the function values (scaled for a maximum of 1. 0) corresponding to the 
breakpoints contained in NRMPN with Y equal to 1. 35. While the input to the sub- 
routine is the normalized pressure ratio, it need not be included as an argument in the 
call statement when following a similar call to CLCFN. The output of the subroutine is 
the scaled-fraction flow function. A FORTRAN listing of CLCWN is given in appendix E. 

NMXTR: Subroutine NMXTR converts a scaled- fraction pressure ratio to a format 
suitable for use with the library- supplied function generation routines. NMXTR is 
called by subroutine TRAT and contains the pressure ratio breakpoints (scaled for a 
maximum of 15. 0) used in the computation of the isentropic temperature rise parameter 
for fans and compressors. The input to the subroutine is the rescaled pressure ratio. 
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The outputs of the subroutine are the normalized pressure ratio and an integer ERRFLG 
used to signal an out- of- range input. A FORTRAN listing of NMXTR is given in appen- 
dix E. 

CLCTR: Subroutine CLCTR generates the isentropic temperature rise parameter 
for a fan or compressor. CLCTR is called by subroutine TRAT and con- 
tains the parameter values (scaled for a maximum of 1. 25) corresponding to the break- 
points contained in NMXTR with Y equal to 1. 35. The input to the subroutine is the 
normalized pressure ratio from NMXTR. The output of the subroutine is the scaled- 
fraction temperature rise parameter. A FORTRAN listing of CLCTR is given in appen- 
dix E. 

MAPFUN: Subroutine MAPFUN generates functions of two variables by radial 
interpolation. MAPFUN is called by subroutines FNCMP and TURB to calculate perfor- 
mance map outputs for specified pressure ratios and rotor speeds. The inputs to 
MAPFUN are a map index MN, the scaled-fraction pressure ratio, and the scaled- 
fraction corrected speed parameter. The map index is defined as follows: MN=1 for a 
fan-hub flow map, MN=2 for a fan-hub efficiency map, MN=3 for a high-pressure- 
compressor flow map, MN=4 for a high-pressure-compressor efficiency map, MN=5 for 
a high-pressure-turbine flow map, MN=6 for a high-pressure-turbine enthalpy drop map, 
MN=7 for a low-pressure-turbine flow map, MN=8 for a low-pressure-turbine enthalpy 
drop map, MN=9 for a fan-tip flow map, and MN=10 for a fan-tip efficiency map. The 
output of the subroutine is the appropriate scaled map variable. The scaled map data 
and search indices are shared by MAPFUN through the use of COMMON arrays. A 
FORTRAN listing of MAPFUN is given in appendix E. 

MOOR: Subroutine MOOR signals the occurrence of out-of-range inputs to maps 
generated by MAPFUN, Termination of input data for fans and compressors at the stall 
line results in MOOR signaling fan or compressor stall. Subroutine MOOR may be 
written to suit the needs of the user. A FORTRAN listing of the existing version of 
MOOR in given in appendix E. 

In addition to the subroutines previously described, the digital program makes use of 
a number of library- supplied routines. Subroutines VBNS, FNGNl, FNLKl, and CBNS 
are part of the library function generation system (ref. 11) and are used to generate 
functions of one variable. These routines are called by subroutines CLCFN, CLCWN, 
and CLCTR. For other computing systems, similar function generation routines are 
usually available. Library routines are also required for performing scaled-fraction 
arithmetic, integer and logical operations, input-output control, and hybrid linkage. 
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Analog Program 


The analog portion of the HYDES program is used, primarily, to perform the inte- 
gration with respect to time associated with the engine dynamics. Appendix F is a 
summary of the scaled equations, solved on the analog computer. For each intercom- 
ponent volume, the scaled- fraction stored mass, temperature, and pressure are compu- 
ted by using the scaled versions of equations (1) to (3). Figure 11 illustrates the analog 
calculation of pressure P„ and temperature T„. Other pressures and temperatures 
are calculated in the same way. 

The analog coefficients C(i) are calculated by the user from the definitions contained 
in appendix D. Scaled-fraction rotor speeds are computed by using the scaled version 
of equation (31). 

Inputs to the analog (through the DAC's) are the scaled stored-mass derivatives 
DW , the scaled temperature derivative terms T ^DN, and the scaled torque differen- 

tials DTQTl and DTQT2. Derivatives for volumes or spools not being used are set to 
zero by the digital program (see main program listing). Control lines are also set to 
allow the user to hold the appropriate analog integrators in the initial condition (IC) if 
drift becomes a problem. Also input to the analog (for display purposes) are the high- 
pressure-compressor map pressure ratio PRC and flow parameter FPC, total engine 
thrust TAUT, and either the fan-hub flow WDFl for HBPR=’’FALSE" or the fan-tip flow 

ADC ^ 


+ P3 



Figure 11. - Analog computation of high -compressor discharge conditions in 
combustor inlet volume. 
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WDF2 for HBPR="TRUE. " The choice of variables to be transmitted to the analog for 
display was dictated by the shortage of available DAC channels and might be changed by 
the user for a particular problem. The outputs from the analog (through the DAC's) are 

the scaled pressures P , temperatures T , and the fraction of design rotor speeds 

PCNTl and PCNT2. Also supplied to the digital by the analog are the scaled nozzle 
areas ANl and AN2, the scaled inlet pressure P2, and fuel flow WDF . Again, the 
choice of these variables was dictated by the shortage of available ADC channels. 

The basic HYDES program does not provide for the simulation of control system 
dynamics. The shortage of hybrid interface channels does, however, dictate that the 
simulation of controls take place on the analog computer. The simulation of a selected 
engine and fuel control is described in a following section. That example illustrates 
the relation between an analog- simulated fuel control and the basic hybrid engine simula- 
tion and demonstrates the transient operation of the program. The analog portion of the 
simulation (without controls) requires 16 summers, 20 integrators, 18 quarter -square 
multipliers, and 57 attenuators. 


Input Data Preparation 

The use of the HYDES program requires that the user provide input data appropriate 
to the engine system being simulated. These data include (1) values for the logical 
variables HBPR, MIX, STRM3, TURBJl, TURBJ2, SUPER, and SPLIT which define 
the engine configuration; (2) values for the digital coefficients SC(i); (3) initial condi- 
tions AI(i) for the DAC variables; (4) unsealed performance data for the fans, compres- 
sors, and turbines; (5) scale factors to be used in scaling the component performance 
data; and (6) values for the analog coefficients C(i). The following sections describe 
the preparation of this input data and the steps required to input the data to the computer. 

Configuration specification . - The execution of the data input program and the main 
digital program (appendix E) requires that the user specify the engine configuration to 
be simulated. For the existing programs, the configuration specification is accomplished 
by depressing sense switches at the control console. Table IV summarizes the configu- 
ration specification procedure. 

Prior to executing the data input program, the user must position sense switches 2, 

5, 6, and 7 as dictated by table IV. For the one-spool turbojet case (configuration G), 
depressing sense switch 5 results in the logical variable TURBJl being set "TRUE" in 
the data input program. By sensing this value of TURBJl, the data input program 
accepts performance data for only one compressor and turbine and stores the resultant 
scaled component data in the proper locations. For configurations A, B, E, and F, 
sense switch 2 is depressed by the user, resulting in the logical variable HBPR being 
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TABLE IV. - CONFIGURATION SPECIFICATION PROCEDURE FOR THE EXISTING HYDES PROGRAM^ 
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^ denotes ’’yes’’ and N denotes ”no. 
^SSW denotes sense switch. 























set "TRUE" in the data input program. This provides for the input and storage of per- 
formance data for a fan or compressor in the bypass stream (fig. 5). For either configu- 
ration E or F, the user may depress sense switch 7, resulting in the variable SUPER 
being set "TRUE" in the data input program. This option (-2) bypasses the data input 
and storage required for mapping the fan hub (low-pressure compressor). A linear 
relation between the "supercharger" pressure ratio and corrected speed is used in the 
main program in place of the usual performance map. 

The execution of the data input program results in the required component perfor- 
mance data, digital coefficients, and DAC initial conditions being stored in the proper 
storage locations for the selected configuration. However, further configuration defini- 
tion is required prior to executing the main program (table IV). If mixing of streams is 
desired (configurations B, D, or F), sense switch 3 is depressed, resulting in a logical 
variable MIX being set "TRUE” in the main program. This variable is tested in the 
program to accomplish the necessary branching and computation associated with volume 
V,^ and the mixed nozzle. For the three-stream configurations (configurations A or B), 
sense switch 4 is depressed, resulting in the logical variable STRM3 being set "TRUE. " 
This option provides for the branching and computation associated with volume ^2 2 
and the fourth nozzle (fig. 1). For configurations A, B, E, and F, depressing sense 
switch 8 results in the variable SPLIT being set "TRUE" in the main program. The 
program then interprets the fan-tip data as being in terms of the total fan flow (-1 
option). For this case, the flow output of the N=3 call to subroutine FNCMP (fan tip) 
is assumed to be the total fan flow. The bypass flow is then calculated by subtracting 
the fan-hub flow from the total fan flow. Since the torque output of the N=3 call to 
FNCMP is computed by using the total fan flow, it must be multiplied by the ratio of 
bypass flow to total fan flow. 

Digital coefficients and DAC initial conditions . - The solution of the scaled equations 
in the main program (appendix E) requires that the digital coefficients SC(i) be stored in 
a COMMON array by using the data input program. The data input program is set up to 
accept NSC values of SC(i), where NSC < 125. The existing program uses 103 values 
of SC(i). The addition of a limited number of new equations and coefficients to the main 
program may be accomplished without modifying the data input program. The values of 
the SC(i)'s are calculated by the user from the definitions given in appendix D. The 
ordering of calculations in the program was chosen to allow the SC(i)'s to be represented 
as scaled fractions. If a calculated SC(i) equals or exceeds unity, the coefficient must 
be redefined or scale factors changed to correct the situation. 

In addition to the NSC values of the digital coefficients, the user must also specify 
initial conditions to be transmitted to the analog through the 24 DAC’s. Since, for the 
existing program, the DAC's are used only for transmitting derivatives and display 
variables, the initial conditions may be set to zero. 
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In the Lewis hybrid computing system, digital input is accomplished by using paper 
tape and a high-speed paper-tape reader. Table V illustrates the format used for the 
input of the digital coefficients and DAC initial condition data. The first line (starting 
in column 1) contains the integer value of NSC (14 format). The next lines contain the 
scaled-fraction values of the digital coefficients (5S8 format). The last five lines con- 
tain the scaled- fraction values of the DAC initial conditions (5S8 format). 

For certain engine configurations, a particular digital coefficient may not be needed. 
A value of unity (0. 99999) should be inserted in the data. 

Component performance data. - The data input program is used to input, scale, and 
store performance data for the fans, compressors, and turbines associated with the 
engine being simulated. Fan and compressor data, in the form shown in figure 9, are 
required. Turbine data must be in the form shown in figure 10. 

Fan and compressor performance data are usually found in the form shown in figure 
12. The conversion of this data to the desired form merely involves reversing the 
pressure ratio and corrected flow axes for the flow map and the reading of efficiency 
values along constant speed lines at selected values of pressure ratio for the efficiency 
map. 


TABLE V. - DIGITAL COEFFiaENT AND DAC INITIAL 
CONDITION INPUT DATA (CONFIGURATION E) 
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Integer value of NSC 

.51867 

. 00000 

. 00000 

.99999 

.99999 



.45731 

. 59066 

.66125 

. 40000. 

, 99999 



.31877 

. 40000 

.06667 

.03333 

. 80000 



.38623 

. 80000 

.43427 

. 16491 

. 26667 



. 00000 

. 00000 

. 00000 

. 10000 

.99999 



. 00000 

.07500 

. 00000 

. 00000 

. 00000 



. 00000 

. 00000 

. 00000 

. 00000 

. 00000 



. 00000 

. 13906 

. 20000 

. 06667 

.03333 


^ Scaled-fraction values of 

. 02000 

.99158 

. 14615 

. 02000 

.03333 


digital coefficients 

,02667 

.99999 

. 00000 

. 00000 

. 20000 



. 99999 

.99999 

. 20000 

. 25000 

.40472 



.39203 

.23858 

,37063 

, 00000 

. 00000 



. 12752 

. 17090 

. 10000 

. 12500 

, 29280 



. 40436 

.76349 

.49547 

.63501 

. 96400 



.79212 

.96400 

.62554 

.82422 

.75317 



. 99999 

.54454 

. 99999 

.16667 

.66667 



. 08333 

.04122 

.28516 

,01532 

. 00000 



, 00000 

. 00000 

. 22443 

. 00000 

. 00000 



. 00000 

-.10539 

. 00045 

- . 00801 

. 25409 



. 24941 

. 25028 

. 50000 

. 28284 

. 00500 

> 


,06667 

. 00000 

. 73500 




. 00000 

. 00000 

. 00000 

. 00000 
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Figure 12. - Typical fan or compressor performance data. 



Figure 13. - Typical turbine performance data. 


Turbine performance data are usually found in the form shown in figure 13. A 
turbine work function dht = Ah/T.^ is plotted against the turbin e cor rected speed for 
constant values of the turbine corrected flow function tff = w V^/Pin and efficiency. 
Since the hybrid program uses pressure ratio as an independent variable in turbine 
calculations, the performance data shown in figure 13 must be converted to a more 
suitable form. 

A digital program, written in FORTRAN for running on the IBM 7094 computer, has 
been developed by the authors of reference 8 to accomplish the turbine map conversion. 
Appendix G contains a listing of that program. For each selected speed parameter (in 
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ascending order), values of the flow function and efficiency are read from figure 13 for 
selected values of the work function. These values, together with the desired pressure 
ratio breakpoints, are read into the conversion program. The conversion program 
proceeds (1) to compute the pressure ratio at each specified point on the map; (2) to 
interpolate, for each speed parameter, the input data to determine values of the work 
function and flow function at the specified pressure-ratio breakpoints; (3) to convert 
the work function and flow function values to the desired flow and enthalpy drop para- 
meters (i. e. , hpt = dht/tnp and fpt = tff/tnp). Since the available turbine map data 
might have to be scaled to the design point of a theoretical engine, the map scaling 
coefficients DHCF and ETTCF must be used in the calculation of actual pressure ratios 
in the conversion program. Also required as input to the conversion program are the 
specific-heat ratio GAM and the specific heat CP for the turbine. The values corres- 
ponding to the design- point discharge temperature and turbine fuel-air ratio are used. 
The output data from the conversion program are plotted in the form shown in figure 10 
to facilitate the removal of any irregularities that may exist in the computed data. 

The configuration specification outlined in table IV enables the user to input only the 
data for components actually required in the simulation. The order of input of map data 
is (1) fan-hub flow map (if required), (2) fan-hub efficiency map (if required), (3) high- 
pressure-compressor flow map, (4) high-pressure-compressor efficiency map, (5) 
high-pressure-turbine flow map, (6) high-pressure-turbine enthalpy drop map, (7) low- 
pressure-turbine flow map (if required), (8) low-pressure-turbine enthalpy drop map 
(if required), (9) fan- tip flow map (if required) (in terms of total fan flow for SPLIT= 
"TRUE" option), and (10) fan-tip efficiency map (if required). 

For the existing data input program, the number of curves per map NCV and the 
number of points per curve NPT have been set to eight. These values may be changed 
by the user, provided the limits of MAPFUN (NPT < 10, NCV < 8) are not exceeded. 

Table VI illustrates the format used for the input of the component performance 
data. For convenience, all data are placed on one data tape. For each map, the first 
line of data (starting in column 1) contains the map scale factors XSC, YSC, and ZSC 
(3F9. 4 format). These scale factors are real numbers and are used to scale the map 
pressure ratio, speed parameter, and map output, respectively. For example, if the 
maximum pressure ratio on the map is 10. 0, a scale factor of XSC=20 might be chosen. 
The second line of data contains the eight values of the unsealed speed parameter for the 
map in ascending order (8F9, 3 format). The next 16 lines contain the unscaled-pressure- 
ratio map output pairs for each specified speed parameter. The data are arranged 
(four pairs to a line) in order of ascending pressure ratio, starting with the data for the 
lowest speed parameter. An 8F9. 4 format is used. (The scale factor on the speed 
parameter for fans and compressors must be 2. 0, and the scale factor on turbine 
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TABLE VI. - Continued. COMPONENT MAP SCALE FACTORS AND DATA (CONFIGURATION E) 
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TABLE VI. - Continued. COMPONENT MAP SCALE FACTORS AND DATA (CONFIGURATION 

(c) High-pressure turbine 
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TABLE VI, - Continued. COMPONENT MAP SCALE FACTORS AND DATA (CONFIGURATION E) 
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Concluded. COMPONENT MAP SCALE FACTORS AND DATA (CONFIGURATION 
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pressure ratio must be 1. 0. ) The inclusion of scale factors as input data allows the 

user to input dimensioned map data in any desired system of units. 

Analog coefficients. - Prior to the execution of the main program, the analog pro- 
gram should be set up“to perform the calculations listed in appendix F. The amplifier 
gains and attenuator settings should be set to achieve the analog coefficients C(i) 
defined in appendix D. The existing program (without controls) uses 57 values of C(i). 

If for a particular engine configuration, a portion of the analog simulation is not 
used(e.g., calculation of W22, T22, and P22 for STRM3=”FALSE” engines), the 
corresponding coefficients should be set by the user to values that will not result in 
analog component overloads. Since the initial condition on the stored mass is used as 
the denominator in the temperature calculation (fig. 11), it should not be set to zero. 


EXAMPLES 

To demonstrate the flexibility of the hybrid computer program, together with its 
steady-state and transient computing capabilities, two engine types are discussed as 
examples. The first - a single-spool, nonafterburning turbojet (configuration G) - 
represents the minimum requirement in terms of component performance mapping 
capability. The second - a high-bypass, two-spool, two-stream turbofan (configura- 
tion E) - requires the full performance mapping capabiUty of the program. 


Single -Spool Turbojet (Configuration G) 

A single-spool turbojet, described in table VII, was operated at NASA's Lewis 
Research Center to determine the effect of inlet distortion on compressor stall limits. 
The HYDES program was used to support those studies, and results from the program 
will serve to demonstrate both the off-design and dynamic capabilities of the program. 

Table VII lists the design-point and engine parameter data for the selected turbojet 
engine. The unity values for the component scaling coefficients (WACCF, DHHPCF, 
etc. ) shown in table VII indicate that the performance data for the compressor and 
turbine were available and not a result of scaling of other component data. The 
compressor temperature interpolation constant was adjusted to give T^ = Tg at e 
design point by using equations (B20) to (B23). Similarly, the temperature interpolation 
constant for the combustor was adjusted to give 

h^j - 1^1 + (f/a)4jh4 - (f/a)4n^(HVF) (44) 
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TABLE VII. - DESIGN- POINT AND ENGINE PARAMETER DATA FOR SELECTED ONE-SPOOL TURBOJET (CONFIGURATION G) 
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TABLE VIIL - SCALE FACTORS FOR SELECTED ONE-SPOOL 
turbojet (CONFIGURATION G) 

[Computer time, V, 100 t.] 


Unsealed variable, x 


2.1 


^3 

W, 

W, 


w 


tl 

"nl 


Wt 


'^bll 

Hr 

Anl 

^nl- ^tot 
tlnp 

fptl 


hptl 

pre 

fpc 

enp 

ceff 


Scaled variable, X 


P21 

P3 

P4 

P5 

P7 

W3 

W4 

W5 

WDC 

WDB 

WDTl 

WDNl 

WDF 

WDBLl 

TRQTI, TRQC 
ANl 

TAUNl, TAUT 

TINP 

FPTl 

DHTl 

HPTl 

PRC 

FPC 

CNP 

CEFF 


Scale factor, SF^ 


13.790 N/cm (20 psia) 
137.90 N/cm^ (200 psia) 
137.90 N/cm^ (200 psia) 
34.474 N/cm^ (50 psia) 

13. 790 N,/cni^ (20 psia) 
0.22680 kg (0. 5 Ibm) 
0.09072 kg (0.2 Ibml 
0.09072 kg (0. 2 Ibm) 
45.359 kg/sec ( 100 Ibm/sec) 


0.45359 kg/sec (1.0 Ibm/sec) 

0. 90718 kg/sec (2. 0 Ibm/ sec) 

4. 5194x10^ N-cm (4xl0‘* in.-lbf) 
1290.3 cm^ (200 in.^) 

2. 2241x10'* N (5000 Ibf) 
1341.6 rpm/K*/^ (1000 rpm/°R*/^) 


0.09137 


N-K-cm 


sec-rpm-N \ sec-rpm-lbf 
2. 1088x10^ j/kg (200 Btu/lbm) 


0.25 


Ibm- R-in. 


3118.7 


('■“ 


Btu 


lbm-°R 


kg-K*/2-rpm 
10 

45. 359 kg/sec (100 Ibm/sec) 
2.0 
1.0 


4/2. 


rpm 


at the design point by using equations (B27) to (B30) ^^ 0 ^ of ^ 

values for /3 and yS. assure a steady-state balance (dXg/dt = 0 and dT^/dt 0) 

energy in voices V„ and when continuity of flow is satisfied. 

Based on the desi|n-point data, scale factors were selected for each variable and 
are summarized in table VIII. As noted previously, certain engine variables have fixe 
scale factors (table II). The specified engine parameters and selected amplitude scale 
factors were used to calculate the values of the 103 digital coefficients SC(i) by means 
of the definitions in appendix D. The previously discussed FORMAT was used to prepare 
the digital coefficient and DAC initial condition data tape. The contents of that tape are 

shown in table IX. 
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TABLE IX, - DIGITAL COEFFICIENT AND DAC INITIAL CONDITION 
INPUT DATA (CONFIGURATION G) 


103 

Integer value of NSC 


, 51867 

, 00000 

. 00000 

. 00000 

. 00000 


. 74996 

. 00000. 

. 00000 

. 00000 

. 00000 


.45180 

. 20000 

. 02000 

, 00000 

. 00000 


. 64372 

. 00000 

. 00000 

. 00000 

. 40000 


, 00000 

. 00000 

. 00000 

. 99999 

. 00000 


.00000 

. 00000 

. 00000 

. 00000 

. 00000 


. 00000 

. 00000 

. 00000 

. 00000 

. 00000 


. 00000 

. 20338 

. 50000 

. 02000 

. 00000 


, 00000 

. 98028 

. 07277 

. 01000 

. 01000 


. 00000 

. 00000 

. 00000 

. 00000 

. 00000 

^ Scaled-fraction values of 

. 00000 

. 00000. 

. 25000 

. 00000 

. 23334 

digital coefficients 

. 00000 

. 31390 

. 00000 

. 00000 

. 00000 


. 08571 

. 00000 

, 08000 

, 00000 

. 67570 


. 00000 

. 00000 

.99999 

. 00000 

. 00000 


.96400 

. 00000 

. 00000 

. 68046 

. 00000 


. 00000 

. 32117 

. 00000 

, 00000 

. 66667 


. 00000 

. 00000 

. 10012 

. 00000 

. 00000 


. 00000 

. 00000 

. 00000 

. 00000 

. 00000 


. 00000 

. 00000 

. 00000 

. 00000 

. 00000 


. 25000 

. 00000 

. 00000. 

. 00000 

. 00000 


. 00000 

. 00000 

.73500 



. 00000 

. 00000 

. 00000 

, 00000 

. 00000 ^ 


. 00000 

. 00000 

. 00000 

. 00000 

. 00000 


. 00000 

. 00000 

. 00000 

. 00000 

. 00000 

Scaled-fraction value of 

. 00000 

. 00000 

. 00000 

. 00000 

. 00000 

initial conditions 

. 00000 

. 00000 

. 00000 

. 00000 

J 



The compressor and turbine performance data, in the form required by the hybrid 
program, are available in reference 8. Table X shows the contents of the prepared 
component performance data tape. 

The volumes listed in table VII resulted in a minimum time constant (eq. (42)) of 
2. 6 milliseconds. A time scale factor SF^’ of 100 was selected to give a minimum ratio 
of digital to analog frequencies equal to 70. 8. The amplitude and time-scale factors 
were used in the calculation of the analog coefficients C(i) as defined in appendix D. 

The resulting coefficients were implemented by using amplifiers and attenuators and are 
listed in table XI. 

Prior to loading and executing the data input program, configuration G was specified 
as per table IV. Sense switch 5 was depressed for the single-spool turbojet case. The 
data input program (appendix E) prompts the user to load (in order) the digital coeffi- 
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TABLE X. -COMPONENT MAP SCALE FACTORS AND DATA (CONFIGURATION G) 
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TABLE X. > Concluded. COMPONENT MAP SCALE FACTORS AND DATA (CONHCURATION G) 
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TABLE XI. - CALCULATED ANALOG 
COEFFICIENTS FOR SELECTED 
TURBOJET ENGINE (CONFIGURATION G) 


■ 


■ 


1 

0 

30 

0 

2 

.9999 

31 

. 9999 

3 

0 

32 

0 

4 

.5187 

33 

0 

5 

.0000 

34 

. 1362 

6 

.9999 

35 

. 9999 

7 

0 

36 

. 3542 

8 

0 

37 

, 0187 

9 

.7472 

38 

. 1969 

10 

.4843 

39 

, 5000 

11 

. 2000 

40 

0 

12 

. 2000 

41 

. 2732 

13 

.5636 

42 

0 

14 

.4112 

43 

0 

15 

.5000 

44 

0 

16 

.5000 

45 

0 

17 

.5851 

46 

0 

18 

.4102 

47 

0 

19 

.5000 

48 

0 

20 

.5000 

49 

. 2075 

21 

.9999 

50 

0 

22 

0 

51 

0 

23 

0 

52 

0 

24 

0 

53 

0 

25 

0 

54 

0 

26 

. 1417 

55 

.7350 

27 

.9999 

56 

.5260 

28 

0 

57 

. 7030 

29 

0 




dents, DAC initial conditions, and component performance data tapes. After the user 
loads the data listed in tables K and X, the data are stored in the appropriate COMMON 
blocks. 

As indicated in table IV, no further steps are required to define the single- spool 
turbojet configuration prior to loading and executing the main digital program. Prior to 
reading ADC values and calculating derivatives, the main digital program outputs initial 
conditions for the derivatives (through the DAC’s) and puts the analog console in the 
initial condition (IC) mode. The design-point values of pressure, temperature, nozzle 


46 




area, and fuel flow rate are computed by using the preset analog coefficients and are 
transmitted to the digital through the 24 ADC's. The stored digital coefficients result 
in the digital program computing near-zero derivatives for the temperatures and stored 
masses, which are transmitted to the analog through the 24 DAC's. This balanced 
condition results in minimal changes in the engine variables when the analog is manually 
placed in the "operate" mode and the integrations commence. 

To allow transient variations in the simulated turbojet's operating point, a simplified 
fuel control model was added to the analog portion of the simulation. A proportional- 
plus-integral speed controller was assumed. The unsealed equation, describing the 
controller, is 


4wj, = - N,j).K2 (^5) 

Gain values of Kj = 2.211x10"^ kg/sec- rpm (9.306x10"^ Ibm/sec-rpm) and 
Kg = 5.620x10"'* kg/sec^-rpm (1.239x10"^ Ibm/sec^-rpm) were selected to give satis- 
factory response characteristics. The scaled version of equation (45) was implemented 
on the analog to provide the scaled fuel flow WDF to the digital program through an ADC. 

w 

With the fuel controller maintaining the compressor speed at the design-point value, 
a steady^state listing of selected engine variables (scaled) was obtained at the teletype 
by depressing sense switch 1 (appendix E). The resulting listing is shown in table XII, 
together with the corresponding values of the unsealed variables. A comparision of 
tables Xn and VII indicates the following maximum differences between the computed 
and specified design points; pressures, 0. 27 percent; temperatures, 0. 30 percent; 
flow rates, 0. 56 percent; thrust, 0. 44 percent. 

For the Lewis turbojet tests, a baseline stall limit was to be established, experi- 
mentally, by reducing the core nozzle area Anl at selected rotor speeds until stall was 
observed. Difficulties arose, however, at high speeds because turbine temperature 
limits were reached before stall occurred. The possibility of decreasing the stall mar- 
gin by decreasing the turbine nozzle area was investigated by using the hybrid simula- 
tion. The experimental tests were to be conducted with a compressor inlet pressure of 

O 

6. 895 N/cm (10 psia). The hybrid- simulated inlet pressure was reduced by decreasing 
the analog coefficient C(26) (appendix D) and thus causing a shift in the 100-percent-speed 
operating point. Table XIII shows the steady- state computer listing for the reduced in- 
let pressure. 

Stall data were obtained for both the nominal turbine area and an area equal to 85. 7 
percent of the nominal area. The decreasing of the turbine nozzle area was simulated 
by increasing the digital coefficient SC(57) (appendix D). The fuel controller was used 
to maintain the rotor speed at desired values, ranging from 100 percent to 87 percent of 
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TABLE XII. - STEADY-STATE DATA FOR SELECTED TURBOJET ENGINE 
(CONFIGURATION G) AT HIGH- PRESSURE -COMPRESSOR INLET 
PRESSURE OF 10. 133 N/cm^ (14.696 psia) 

[Percent of design speed, 100, ] 

(a) Scaled teletype output 


WDFl 
. 00009 

WDC 

.43911 

WDB 

.42703 

WDTl 

,43124 

WDT2 

.44580 

WDNl 

.44580 

WDN2 
. 00000 

WDN3 
. 00000 

WDN4 
. 00000 

WDF2 
. 00000 

P2 

,73498 

P21 

.73437 

P35 

-.00036 

P3 

. 49206 

P4 

.45483 

P5 

,65283 

P6 

.67283 

P7 

.73474 

P22 

.00024 

P32 
. 00000 

T35 
. 00000 

T3 

.48376 

T4 

.41027 

T5 

.40905 

T6 

,65447 

T7 

.20727 

T2 

.51867 

T21 
. 51818 

TAUNl 

.56680 

TAUN2 
. 00000 

TAUN3 
. 00000 

WDF 

.69836 

PCNTl 

.49938 

PCNT2 
. 00000 

PCNFl 

.00000 

PCNF2 

.00000, 

PCNC 

.49938 

ANl 
. 52600 

AN2 
. 00000 

AN3 
. 00000 

AN4 
. 00000 

ABLC 

.00000 

ABLS 
. 00000 

TAUN4 
. 00000 

TAUT 

.56677 


(b) Unsealed data 


. 19.92 (43.91) 
. 19.37 (42,70) 
. 19. 56 (43.12) 
. 20.22 (44.58) 
, 10. 12 (14.68) 
. 67.85 (98.41) 
. 62.72 (90.97) 
. 22. 50 (32.64) 
. 287.9 (518,2) 
. 537. 5 (967. 5) 
. . 1139 (2051) 
. . 908.9 (1636) 
0.3168 (0.6984) 

0.9988 

. 12 606 (2834) 


High-pressure- compressor flow rate, w^, kg/sec (Ibm/sec) 

Combustor flow rate, w^, kg/sec (Ibm/sec) . 

High-pressure-turbine flow rate, w^p kg/sec (Ibm/sec) 

Core nozzle flow rate, w^^, kg/sec (Ibm/sec) 

High-pressure- compressor inlet pressure, P2 ^ N/cm (psia) . . . . 

Combustor inlet pressure, P^, N/cm^ (psia) 

High-pressure- turbine inlet pressure, P^, N/cm (psia) 

Low-pressure-turbine inlet pressure, P^, N/cm"^ (psia) 

High-pressure- compressor inlet temperature, Tg K (^R) 

Combustor inlet temperature, K (^R) 

High- pressure- turbine inlet temperature, T^, K (°R) 

Low-pressure- turbine inlet temperature, T^, K (°R) 

Fuel flow rate, w^, kg/sec (Ibm/sec) 

Ratio of high- pres sure- turbine rotational speed to design, Ntl'^^tljdes 
Total thrust, N (Ibf) 
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TABLE Xni. - STEADY-STATE DATA FOR SELECTED TURBOJET ENGINE 


(CONFIGURATION G) AT HIGH- PRESSURE-COMPRESSOR INLET 
PRESSURE OF 6.895 N/cm^ (10 psia) 


[Percent of design speed, 100. ] 


(a) Scaled teletype output 


WDFl 
. 00009 

WDC 

.29873 

WDB 

.28720 

WDTl 

.29428 

WDT2 

.30361 

WDNl 

,30361 

WDN2 
. 00000 

WDN3 

.00000 

WDN4 

,00000 

WDF2 
. 00000 

P2 

.73498 

P21 

.49938 

P35 

.00024 

P3 

.33068 

P4 

.31384 

P5 

,46362 

P6 

.46362 

P7 

.73425 

P22 
. 00000 

P32 

-.00024 

T35 
. 00000 

T3 

.48291 

T4 

.41918 

T5 

.42114 

T6 

.67382 

T7 

.20727 

T2 

. 51867 

T21 

.51855 

TAUNl 

.30416 

TAUN2 

.00000 

TAUN3 

,00000 

WDF 

.49951 

PCNTl 

.50012 

PCNT2 

.00012 

PCNFl 
. 00000 

PCNF2 

.00000 

PCNC 
. 50012 

ANl 

.52612 

AN2 
, 00000 

AN3 

.00000 

AN4 

.00000 

ABLC 
. 00000 

ABLS 
. 00000 

TAUN4 

.00000 

TAUT 

.30413 


(b) Unsealed data 


High-pressure- compressor flow rate, w^, kg/sec (Ibm/sec) 

Combustor flow rate, kg/sec (Ibm/sec) 

High-pressure-turbine flow rate, w^p kg/sec (Ibm/sec) 

Core nozzle flow rate, w , kg/sec (Ibm/sec) 

n 1 2 

High-pressure- compressor inlet pressure, ?2 N/cm (psia) 

Combustor inlet pressure, P«, N/cm^ (psia) 

2 

High-pressure-turbine inlet pressure, P^, N/cm (psia) 

^ 2 

Low-pressure-turbine inlet pressure, Pg, N/cm (psia) 

High-pressure- compressor inlet temperature, T 2 K (”r) 

Combustor inlet temperature, Tg. K (°R) 

High-pressure- turbine inlet temperature, T^. K ('^R) 

Low-pressure-turbine inlet temperature, Tg. K (°R) 

Fuel flow rate, Wp, kg/sec (Ibm/sec) 

Ratio of high-pressure- compressor rotational speed to design. N^/Ng (jgg 
Total thrust, Fj.^j., N (Ibf) 


. 13.55 (29.87) 
. 13.03 (28.72) 
. 13.35 (29.43) 
. 13.77 (30.36) 
. 6.886 (9.988) 
. 45.60 (66.14) 
. 43.28 (62.77) 
. 15.98 (23.18) 
. 288.1 (518.6) 
. 536.6 (965.8) 
. . 1164 (2096) 
. . 935. 5 (1684) 
0.2266 (0.4995) 

1.000 

. 6765.7 (1521) 
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Figure 14. - Effect of turbine nozzle area on compressor 
stall margin for turbojet engine {configuration G). 


the design speed. At each selected speed, the core nozzle area was reduced by decreas- 
ing the coefficient C(56) on the analog. Compressor stall was detected by monitoring 
the compressor flow map variables PRC and FPC on an x-y plotter. These variables 
are transmitted to the analog through the DAC's. Oscillations along constant-speed lines 
were observed when stall was initiated. Figure 14 shows the results of that study. With 
the nominal turbine area, stall could not be initiated without exceeding the maximum Tg 
or minimum A^^j limits. For the 100- , 94-, and 87-percent rotor speeds, however, the 
reduced turbine area allowed the compressor to stall with temperatures below the Tg 
limit and nozzle areas above the minimum A^j. With this information, the required 
turbine modifications were made and the experimental program continued. 


Turbofan (Configuration E) 

To further demonstrate the capabilities of the HYDES program, a two-spool, two- 
stream turbofan (configuration E) was selected for simulation. Configuration E was 
selected because of its extensive function generation requirements. A total of 10 func- 
tions of two variables were required to describe the low-pressure compressor (boosted 
fan hub), high-pressure compressor, high-pressure turbine, low-pressure turbine, and 
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bypass fan (fan tip). Table XIV lists the steady-state characteristics of the selected 
turbofan engine at 50 percent of its design thrust. The nonunity values of the component 
scaling coefficients indicate that the performance data being used were scaled from other 
engine components. As in the case of the turbojet, temperature interpolation constants 

> Prof P > and 'wnre adjusted to give steady-state energy balances in volumes 
^2. P V3, 2’ V3, and V^, respectively. 

’ Scale factors were selected for each variable and are summarized in table XV. The 
specified engine parameters and selected amplitude scale factors were used to calculate 
the 103 digital coefficients SC(i). The contents of the digital coefficient and DAC initial 
condition tape are shown in table V. 

The fan, compressor, and turbine performance data for the selected engine were 
available in the forms shown in figures 12 and 13. The turbine data conversion program 
(appendix G) was used to convert the data for each turbine to the required form. Table 
VI shows the contents of the component performance data tape for the selected turbofan. 

The volumes listed in table XIV resulted in a minimum time constant of 4. 75 milli- 
seconds. A time scale factor SF[ of 100 was selected to give a minimum ratio of digital 
to analog frequencies equal to 68. 3. The amplitude and time scale factors were used to 
calculate the analog coefficients C(i). The resulting coefficients are listed in table XVI. 

Before the data input program was loaded and executed, configuration E was speci- 
fied as per table IV. Sense switch 2 was depressed to denote the HBPR="TRUE" case. 
Sense switch 7 was not depressed since the fan hub or low-pressure-compressor had to 
be represented by a full set of data (SU PE R=" FALSE"). After the digital coefficient, 

DAC initial condition, and component performance data were loaded, no further steps 
were required to define the engine. The main program was then loaded. 

To allow transient variations in the simulated turbofan's operating point, a portion of 
a proposed fuel control system was implemented on the analog computer. The features 
of the fuel control system and its implementation are described in appendix H. Basically, 
the control provides for closed- loop control of fan speed under quasi-steady-state 
conditions. 

With the fuel control maintaining the fan speed at the initial value (73. 2 percent of 
design), a steady-state listing of selected engine variables was obtained at the teletype 
by depressing sense switch 1. The resulting listing is shown in table XVII together with 
the corresponding unsealed values. 

A transient was initiated by ramping the demanded low rotor speed N^2, dem 
73. 2 percent value to 100 percent of the design-point value in 50 milliseconds (5 sec of 
computer time). Figure 15 shows the resultant responses of thrust and turbine inlet 
temperature T^, respectively. These plots were obtained by monitoring the outputs of 
the corresponding analog amplifiers with an x-y plotter. An integrator was used to time 
the transient. With the nominal acceleration schedule of fuel flow (eq. (H12)), 63 per- 
cent of the desired thrust change was achieved in about 1. 1 seconds. However, the 
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TABLE XV. - SCALE FACTORS FOR SELECTED TURBOFAN (CONFIGURATION E) 

[Computer time, t’, lOOt.] 
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TABLE XVI. - CALCULATED ANALOG 
COEFFICIENTS FOR SELECTED 
TURBOFAN ENGINE 


(CONFIGURATION E) 


i 

Coefficient, 

i 

Coefficient. 


c, 


Ci 

1 

0.7500 

30 

0 

2 

. 6390 

31 

. 6954 

3 

, 7500 

32 

, 1754 

4 

.5800 

33 

0 

5 

0 

34 

. 1672 

6 

. 9999 

35 

. 9999 

7 

0 

36 

. 3542 

8 

0 

37 

.0042 

9 

.3752 

38 

. 2014 

10 

.5758 

39 

.4438 

11 

.0750 

40 

, 0139 

12 

.0750 

41 

. 1802 

13 

. 3849 

42 

. 3660 

14 

.4446 

43 

. 1319 

15 

. 1000 

44 

. 2000 

16 

. 1000 

45 

.5387 

17 

. 6284 

46 

. 2000 

18 

.4195 

47 

. 5392 

19 

. 3000 

48 

0 

20 

. 3000 

49 

. 2075 

21 

.8940 

50 

0 

22 

.5206 

' 51 

0 

23 

. 1500 

52 

0 

24 

. 1500 

53 

0 

25 

. 2000 

54 

. 6684 

26 

. 1087 

55 

. 7238 

27 

. 6810 

56 

. 7358 

28 

. 1754 

57 

.3358 

29 

. 2000 






TABLE XVII. - STEADY-STATE DATA FOR SELECTED TURBOFAN ENGINE 
(CONFIGURATION E) AT 73.2 PERCENT OF DESIGN SPEED 
(a) Scaled teletype output 


WDFl 

.43469 

WDC 

.43475 

WDB 

.39212 

WDTl 
. 40072 

WDT2 

.42034 

WDNl 

,44061 

WDN2 
. 52441 

WDN3 

.00000 

WDN4 
. 00000 

WDF2 
. 52346 

P2 

.72378 

P21 

,40197 

P35 

,64294 

P3 

.36254 

P4 

.34692 

P5 

.47534 

P6 

,61486 

P7 

.73474 

P22 

-.00048 

P32 

.65649 

T35 

.53857 

T3 

, 57617 

T4 

.44458 

T5 

.41662 

T6 

. 51660 

T7 

,20727 

T2 

. 51867 

T21 

.57958 

TAUNl 

.16723 

TAUN2 

.24221 

TAUN3 

.00000 

WDF 

.33044 

PCNTl 

.44409 

PCNT2 

.36584 

PCNFl 

.36581 

PCNF2 

.36581 

PCNC 

.44409 

ANl 

.73547 

AN2 

.66809 

AN3 

.00000 

AN4 
. 00000 

ABLC 
. 00000 

ABLS 
. 00000 

TAUN4 
. 00000 

TAUT 

.25891 


(b) Unsealed data 

Low-pressure- compressor flow rate, Wq, kg/sec (Ibm/sec) 29.57 (65.20) 

High-pressure- compressor flow rate, w^, kg/sec (Ibm/sec) 29.57 (65.20) 

Combustor flow rate, kg/sec (Ibm/sec) 26.68 (58.82) 

High-pressure- turbine flow rate, w^^, kg/sec (ibm/sec) 27.26 (60.10) 

Low-pressure- turbine flow rate, Wj. 2 , kg/sec (ibm/sec) 28.59 (63.04) 

Core nozzle How rate, w^^^, kg/sec (ibm/sec) 29.98 (66.09) 

Fan flow rate, Wj 2 ‘. third-stream nozzle flow rate, w^ 2 ’ kg/sec (Ibm/sec) . . . . 475.8 (1049) 

Fan inlet pressure, P 2 , N/cm^ (psia) ^ 9.984 (14.48) 

High-pressure- compressor inlet pressure, P 2 N/cm (psia) 13.86 (20.10) 

Combustor inlet pressure, P 3 , N/cm^(psia) . . 124.9(181.2) 

High-pressure- turbine inlet pressure, P,j, N/cm 119.56 (173.4) 

Low-pressure-turbine inlet pressure, Pg, N/cm 32.77 (47.53) 

Core duct pressure, Pg, N/cm^(psia). • ■ 10.60 (15.37) 

Bypass fan discharge pressure, Pg g-N/cm 11.31(16.41) 

Bypass fan duct pressure, P, g, N/cm^ (psia) 11.08 (16.07) 

Fan inlet temperature, T 2 , k'(*R) ^ 288.2 (518.7) 

High-pressure- compressor inlet temperature, T 2 K (°R) 322.0 (579.6) 

Combustor inlet temperature, Tg, K (°R) 640.0 (1152) 

High-pressure- turbine inlet temperature, T^, K (°R) 1235 (2223) 

Low-pressure- turbine inlet temperature, Tg, K (°R) 925. 5 (1666) 

Core duct temperature, Tg, K (°R) ® (1292) 

Fuel How rate, Wp, kg/sec (Ibm/sec) 0.4496 (0.9912) 

Ratio of high-pressure- turbine rotational speed to design, Ntl/Ntl, des ■ ’ 6882 

Ratio of low-pressure- turbine rotational speed to design, N^2/N^2 des 

Total thrust, F^^^, N (Ibf) • • • • S'? 582 (12 945) 
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High-pressure turbine inlet temperature, 



(a) Thrust. 



(b) Temperature. 

Figure 15. - Response of selected turbofan engine to 50-millisecond ramp demand in low rotor speed. 
Initial speed demand, 73. 2 percent of design; final speed demand, 100 percent of design. Analoq- 
simulated fuel control. Variable coefficient in acceleration schedule 






TABLE XVITI. - STEADY-STATE DATA FOR SELECTED TURBOFAN ENGINE 
(CONFIGURATION E) AT 100 PERCENT OF DESIGN SPEED 
(a) Scaled teletype output 


WDFl 

.63925 

WDC 

.63973 

WDB 

.57156 

WDTl 
. 58990 

WDT2 

.62103 

WDNl 

.64282 

WDN2 

.72888 

WDN3 
. 00000 

WDN4 
. 00000 

WDF2 

.72671 

P2 

.72378 

P21 
. 51513 

P35 

.69653 

P3 

,61328 

P4 

. 58007 

P5 

.77160 

P 6 

.64953 

P7 

.73474 

P22 

-.00036 

P32 

.72326 

T35 
. 55590 

T3 

.68188 

T4 

.57385 

T5 

. 54235 

T 6 

.62341 

T7 

.20727 

T2 

. 51867 

T21 

.63146 

TAUNl 

.43151 

TAUN2 

.47259 

TAUN3 
. 00000 

WDF 
. 72692 

PCNTl 
. 50024 

PCNT2 

.49951 

PCNFl 

.49948 

PCNF2 

.49948 

PCNC 
. 50024 

ANl 

.73559 

AN2 

.66809 

AN3 
. 00000 

AN4 
. 00000 

ABLC 
. 00000 

ABLS 
. 00000 

TAUN4 
. 00000 

TAUT 

.51571 


(b) Unsealed data 


Low-pressure-compressor flow rate, kg/sec (Ibm/sec) 43.51 (95.92) 

High-pressure- compressor flow rate, w^, kg/sec (Ibm/sec) 43.51 (95.92) 

Combustor flow rate, Wj^, kg/sec (Ibm/sec) 

High-pressure- turbine flow rate, kg/sec (Ibm/sec) 40.13 

Low-pressure- turbine flow rate, Wj.2, kg/sec (Ibm/sec) 42.25 (93.15 

Core nozzle flow rate, w^^, kg/sec (Ibm/sec) . • • 

Fan flow rate, Wf«; third-stream nozzle flow rate, w^2^ kg/sec (Ibm/sec) 660. 

2 f 9.984 (14,48) 

Fan inlet pressure, P 2 . N/cm (psia) 17 7 r 

High-pressure- compressor inlet pressure, P 2 N/cm (psia) 

Combustor inlet pressure, P3, N/cm^ (psia) . • ‘ ^ 

High-pressure- turbine inlet pressure, P^, N/cm (psia) 19 . 

Low-pressure-turbine inlet pressure. Pc, N/cm (psia) ^ 

Core duct pressure, Pg, N/cm^ (psia) . . . 12 46 (18 08) 

Bypass fan discharge pressure, Pg g. N/cm (psia) ^ 

Bypass tan duct pressure, Ps. s. fcm (psia) 288.2 (518.7) 

Fan inlet temperature, T 2 , K ( R) 

High-pressure- compressor inlet temperature, T 2 K ( R) 

Combustor inlet temperature, T 3 , K ( R) . . . . 

High-pressure- turbine inlet temperature, T^, K(^R) * 

Low-pressure-turbine inlet temperature, T^, K ( 

core duct temperatoe. Tj KJ R)^ 9.93 (2. 181) 

Fuel flow rate, w^., kg/sec (Ibm/sec) 

Ratio of high-pressure-turbine rotational speed to design, 0 9990 

Ratio of low-pressure-turbine rotational speed to design, Nj^2^^t2,des 5 

Total thrust, N (Ibf) , 
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transient was terminated automatically after 1. 2 seconds when the boosted- fan-hub map 
inputs went out of range (booster stall). The stall condition was caused by the excessive 
turbine inlet temperature shown in figure 15(b). 

To demonstrate the effects of reducing the acceleration schedule of fuel flow, the 
coefficient (eq. (H12)) was reduced in steps of 1.827x10"'^ kg-cm^sec-N 
(1. 0 Ibm/hr-psia). The engine was accelerated successfully with the acceleration 
schedule reduced by more than 3. 654x10"^ kg-cm^sec-N (2. 0 Ibm/hr-psia). Figure 15 
shows the trade-off between reduced turbine inlet temperatures and fast thrust response. 

At the end of one of the transients, a teletype listing of steady-state variables was 
obtained by depressing sense switch 1. The resulting listing is shown in table XVin. 

A comparison of table XVIII with steady-state data generated with GENENG II (ref. 4) 
at this operating point indicates the following maximum differences between the two 
programs: pressures, 0.07 percent; temperatures, 0. 42 percent; flows rates, 0.77 
percent; thrust, 4. 4 percent. 


CONCLUDING REMARKS 

This report has described HYDES, a hybrid computer program capable of simulating 
either one-spool turbojet, two-spool turbojet, or two-spool turbofan engine dynamics. 
The program is also capable of simulating two-or three- stream turbofans with or with- 
out mixing of the exhaust streams. HYDES was developed for running on the Lewis 
Research Center’s Electronic Associates (EAI) 690 Hybrid Computing System. The 
hybrid computer combines the precision and logic capabilities of the digital computer 
with the integration and data output capabilities of the analog computer. The HYDES 
program is intended to eliminate the need for developing individual simulations for 
each new engine study. 

In the HYDES program, the analog computer is used, primarily, for performing 
integration with respect to time. The use of the digital computer to perform all function 
generation and most of the algebraic calculations results in digital cycle times between 
23 and 44 milliseconds. These long cycle times make real-time simulation impossible 
and require time scaling of the dynamics on the analog computer. 

The documentation of the HYDES program, together with the digital computer soft- 
ware (available from the author upon request), should allow the user to quickly imple- 
ment the program on a similar computer. It is expected that this report will also prove 
valuable in the development of both generalized and specific engine simulations on other 
computers. The scaling techniques, discussed in the report, are applicable to all 
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analog simulations. The branching techniques and program structure associated with 
the generalization of the HYDES program should be applicable to all digital simulations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 25, 1973, 

501-24. 
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APPENDIX A 


UNSCALED-VARIABLE SYMBOLS 
2 2 

A effective cross-sectional area, cm (in. ) 

a supercharger map intercept 

b supercharger map slope 

nozzle velocity coefficient 

CNCF general turbine speed scaling coefficient 

CNHPCF high-pressure-turbine speed scaling coefficient 

CNLPCF low-pressure -turbine speed scaling coefficient 

c specific heat at constant pressure, J/kg-K (Btu/lbm-°R) 

specific heat at constant volume, J/kg-K (Btu/lbm-°R) 

ceff high-pressure- compressor map efficiency 

cnp high-pressure-compressor map speed parameter 

ctr high-pressure- compressor temperature rise parameter 

DHCF general turbine enthalpy scaling coefficient 
DHHPCF high-pressure-turbine enthalpy scaling coefficient 

DHLPCF low-pressure-turbine enthalpy scaling coefficient 

dht general turbine map work function, J/kg-K (Btu/lbm-°R) 

ETACCF high-pressure-compressor efficiency scaling coefficient 

ETACF general fan or compressor efficiency scaling coefficient 

ETAFCF fan efficiency scaling coefficient 

ETAICF low-pressure-compressor efficiency scaling coefficient 
ETTCF general turbine efficiency scaling coefficient 

eff general component map efficiency 

F thrust, N(lbf) 

AF thrust increment, N(lbf) 

Fi nozzle function, i = 1 for flow, i = 2 for thrust 

f- functional relation 
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f/a 

fcnp 

fleff 

flnp 

fltr 

fpc 

fpfc 

fpfl 

fpf2 

fpt 

fptl 

fpt2 

f2eff 

f2np 

f2tr 

HVF 

h 

Ah 

hpt 

hptl 

hpt2 

I 

J 


local fuel-air ratio 

general fan or compressor map speed parameter 

low-pressure-compressor map efficiency 

low-pressure-compressor map speed parameter 

low-pressure-compressor ideal temperature rise parameter 

high-pressure-compressor map flow parameter, kg/sec (Ibm/sec) 

general fan or compressor map flow parameter, kg/sec (Ibm/sec) 

low-pressure-compressor map flow parameter, kg/sec (Ibm/sec) 

fan map flow parameter, kg/sec (Ibm/sec) 

general turbine map flow parameter, 

kg-K-cm^/N-rpm-sec (lbm-°R-in. ^Haf-rpm-sec) 

high-pressure-turbine map flow parameter, 

kg-K-cm^/N-rpm-sec (lbm-°R-in. VlM-rpm-sec) 

low-pressure-turbine map flow parameter, 

kg-K-cm^/N-rpm-sec (Ibm-^R-in. Vlbf-rpm-sec) 

fan map efficiency 

fan map speed parameter 

fan map ideal temperature rise parameter 

2 2 \ 
gravitational conversion factor, 100 cm-kg/N-sec (386.3 Ibm-in./lbf-sec ) 

heating value of fuel, 4. 302x10^^ J/kg (18 500 Btu/lbm) 

enthalpy, J/kg, (Btu/lbm) 

general turbine enthalpy drop, J/kg (Btu/lbm 

general turbine map enthalpy drop parameter , 
j/kg-K^^^-rpm (Btu/lbm-°R^^^- rpm) 

high-pressure- turbine map enthalpy drop parameter, 
j/kg-K^"^^-rpm (Btu/lbm- °R^^^-rpm) 

low-pressure-turbine map enthalpy drop parameter, 

J/kg-K^/^-rpm (Btu/lbm-°R^/^-rpm) 

2 2 

polar moment of inertia, N-cm-sec (in.-lbf-sec ) 
mechanical equivalent of heat, 100 N-cm/J (9339. 1 in. -Ibf/Btu) 
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*^f2 

K. 

1 

L 

AL 

MM 

N 

NN 

P 

P 

s 

PRCCF 

PRCF 

PRFCF 

PRICF 

pr 

prc 

prfl 

prf2 

prtl 

prt2 

R 

dt 

S 

T 

AT 


low-pressure-compressor/low-pressure-turbine speed ratio 
fan/low-pressure-turbine speed ratio 
controller gain 

acceleration schedule coefficient, kg-cm /sec-N (Ibm-hr-psia) 

torque, N-cm (in.-lbf) 

differential torque, N-cm (Ibf-in.) 

number of streams fed by control volume 

rotational speed, rpm 

number of streams feeding control volume 

2 2 

total pressure, N/cm (Ibf/in. ) 
static pressure, N/cm (Ibf/in. ) 

high-pressure-compressor pressure- ratio scaling coefficient 
general fan or compressor pressure-ratio scaling coefficient 
fan pressure-ratio scaling coefficient 

low-pressure- compressor pressure- ratio scaling coefficient 
map pressure ratio 

high-pressure-compressor map pressure ratio 

low-pressure-compressor map pressure ratio 

fan map pressure ratio 

high-pressure-turbine pressure ratio 

low-pressure-turbine pressure ratio 

gas constant, N-cm/kg-K (in. -lbf/lbm-°R) 

gas constant of air, 44.83 N-cm/kg-K 
(640.1 in. -Ibf/lbm-^R) 

2 2 2 

pressure loss coefficient, N-sec /kg -cm 

Laplacian operator, sec'^ 

total temperature, K(°R) 

fan or compressor temperature rise, K(°R) 



TFCF 

TFHPCF 

TFLPCF 

t 

dt 

tff 


tnp 

tr 

tlnp 


t2np 

V 
W 

WACCF 

WACF 

WAFCF 

WAICF 

w 

AWj 

X 

y 

a 

P 

V 
e 

V 

T 

<P 


general turbine flow scaling coefficient 
high-pressure-turbine flow scaling coefficient 
low-pressure-turbine flow scaling coefficient 
time, sec 

differential time, sec 

1/2 2 

general turbine map flow function, kg-K ' -cm /N-sec 
(lbm-°R^/^-in. Vlbf-sec) 

general turbine map speed parameter 

general fan or compressor temperature rise parameter 

high-pressure-turbine map speed parameter, 
rpm/K^^^ (rpm/°R^^*'^) 

low-pressure-turbine map speed parameter, 

rpm/K^'^^ (rpm/^R^'^^) 

3 3 

volume, cm (in. ) 

stored mass, kg (Ibm) 

high-pressure- compressor flow scaling coefficient 
general fan or compressor flow scaling coefficient 
fan flow scaling coefficient 

low-pressure-compressor flow scaling coefficient 

flow rate, kg/sec (Ibm/sec) 

output of fuel controller, kg/sec (Ibm/sec) 

general unsealed variable 

position signal 

power lever angle 

temperature interpolation constant 

specific-heat ratio 

control error signal 

efficiency 

time constant, sec 

acceleration schedule outjnit, kg-cm / sec-N (Ibm/hr-psia) 
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frequency, sec 


fjj 


-1 


Subscripts: 

a analog 

accel acceleration 

b combustor 

blc control bleed 

bis interstream bleed 

bll high-pressure-turbine cooling bleed 

bl2 low-pressure-turbine cooling bleed 

c high-pressure compressor 

cr critical 

d digital 

dem demanded value 

des design value 

F fuel 

FB feedback 

f final value 

fm fan map 

fl low-pressure compressor 

f2 fan 

i initial value 

id ideal 

in input to control volume 

j index on flows into control volume 

k index on flows out of control volume 

m measured 

max maximum value 



min minimum value 

nom nominal value 

nl core nozzle 

n2 third- stream nozzle 

n3 mixing nozzle 

n4 second- stream nozzle 

out out of control volume 

ovb overboard bleed 

sc supercharger 

si sea level, 10.133 N/cm^ (14.696 psia) 

std standard day, 288. 1 K (518. 67° R) 

tot total 

tl high-pressure turbine 

t2 low-pressure turbine 

v2 second-stream valve or duct 

v3 third-stream valve or duct 

2 fan inlet 

2. 1 high-pressure-compressor inlet 

2.2 second- stream duct 

3 combustor inlet 

3.2 bypass fan discharge 

3.5 bypass fan duct 

4 high-pressure-turbine inlet 

5 low-pressure-turbine inlet 

6 core duct 

7 mixing volume 

8 nozzle discharge 
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Superscripts: 

(" ) average value 

I 

( ) inlet to control volume 

(• ) time derivative 
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APPENDIX B 


SUMMARY OF UNSCALED EQUATIONS 



PRICF 


fpfl = fj(prfl, flnp) 


(WAICF)(fpfl) 


Wfi - 




r/H 

r 


fleff = l 2 (prfl, flnp) 

= (ETAICF) (fleff) 


fltr = 


(^) 


^fl “ %^"^fl’ 

Tfi = + (1 - ^fi)T2_ 1 


(Bl) 


(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 
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-2.1. 

(BIO) 


(Bll) 

30J(h2 j - h2)wjj 

Lifl 

ITNfi 

(B12) 

l'2 = '4<-2’“) 

(B13) 


N, 


N 


cnp 


c,des 


2.1 


prc 




(B14) 


(B15) 


PRCCF 
= fg(prc, cnp) 


(WACCF)(ft)c) -111 



ceff = fg^rc, cnp) 

= (ETACCF)(ceff) 


(B16) 

(B17) 

(B18) 

(B19) 
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, 30J(h^ -h2 j)w^ 

"c - .K, 

(B25) 

''2.1 ' 1’®^ 

(B26) 

\ = Vs + (* - ^b)T4 

(B27) 


(B28) 

*F 

(t/a)4 = ^ 
*b 

(B29) 

^ iA3-"^4 

" y 

(B30) 


tlnp 


(CNHPCF)N^j 


(B31) 


(B32) 


fptl = f,y(prtl, tlnp) 


(B33) 


hptl = fg(prtl, tlnp) 


^tl = 


(CNHPCF)(N^^)(P 4 )(ft)tl) 

(TFHPCF)(T4) 


(DHHPCF)(hptl)(CNHPCF)(N..) WT. 

Ah = ■ - 

1000 


30J4hj.jW^I^ 






[T 4 , (f/a)J 


[(f/a). 


+1 M 


^bll = "^bll 


Sc ^3^1 As 


is fit by ( 0 . 44208 + 0. 17337 y^). 
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’'’s " 


(CNLPCF)N^2 
t2np = — 

Pfi 

prt2 = — 


fpt2 = fgft)rt2, t2np) 

hpt2 = fjQ(prt2, t2np) 


(CNLPCF)(N^2)(P5)(fpt2) 

(TFLPCF)(T5) 


Ah^2 - 


(DHLPCF)(hpt2)(CNLPCF)(N^2)V^ 

1000 


^2 “ 


30JAhj.2W|.2 


irN, 


t2 


^6 ' ^5 " ^^t2 
hg = £4 [T5, (f/a),] 


(B42) 

(B43) 

(B44) 

(B45) 

(B46) 

(B47) 

(B48) 

(B49) 

(B50) 

(B51) 

(B52) 
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(f/a)0 = 


1 +[ (f/a)5 + l] ^^bl2 


(B53) 


w, io =A) 


''hl2 - bl2 


1^' 


where ^^(Pg/Pg, Vg) is fit by (0.44208 + 0.17337 yg) 


(B54) 


f2np = 


y2, des 


2, des 


(^) . . 


PRFCF 


fpf2 = fjL^(prf2, f2np) 


(WAFCF)(fpf2) 


where for SPLIT="TRUE’', w^g = - w^j; and for SPLIT=”FALSE”, w^g ' ^fm’ 


f2eff = fjg(prf2, f2np) 
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where 


= (ETAFCF)(f2eff) 


f2tr = 




^ f2 " ^3^^f2 ' 


■^f2 " Pf2'^2 + ■ '^f2^ "^3.2 


/ £2tr \ 

T' = I — + 1 1 Tc 


3.2 


^£2 


^ 3 . 2 = ¥^ 3 . 2 * 0 ) 


^£2 “ 


30J(h^^2 - ^2)Wf2 


ttN , 


£2 


/gc ^3 

%vb ~ -^ovb^/ ~ — (~;;~ » "^3 


*^ aV Tg V^a 


>3, Vg) is £it by (0. 44208 + 0,17337 r^) 


w =A / ^2. 1 -^1 / ^8 , Y 

blc ^/T„ , { P, 




A t/ 2.1 \ ^2. 1 


(B60) 

(B61) 

(B62) 

(B63) 

(B64) 

(B65) 

(B66) 

(B67) 


(B68) 
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(B69) 



where is equal to Pg for STRM3="TRUE" and to P,^ otherwise. 

^3.5=f3(T3.5»0) 

'^2. 2 " V^2. 2’ 


(B70) 

(B71) 

(B72) 

(B73) 

(B74) 

(B75) 

(B76) 

(B77) 
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(B78) 




(f/a). = 


(f/a), 


1 + 


[{f/a)( 


+ 1 



where is equal to for HBPR='’TRUE” and to otherwise. 




T. 


‘'nl=C,,nI*nl>gJ'-p, 6'6 '2 ^ p 


s, nl 


'^6 ^ Anl(Ps,nl- P' 


where is equal to 0. 53685 Pg if P 7 /P 0 ^0. 53685 and to P^ if P 7 /P 0 

^„(pr, y ) = 


^n2 ^v, n2^n2 V ^P> 3. 5 3.5 2 \ Pg 1 


/ 


2J - To;. 1-^^^ » '^3.5 )+A„o(P 


n 2 ' s,n 2 


where P^ equal to 0. 53685 Pg g if P 7 /P 3 5 ^ 0 * 53685 and to P,^ if 

P 7 /P 3 g>0. 53685. 


F_o = C„ „o^ 


/ 2J 

V^c 


c„ oT,, 


n3 ~ ^v, n3 ^^n3 W Pj 2 


s, n3 


’ ^7) +^n 3 ^^s,n 3 " ^ 8 ^ 


where P^ ^^3 is equal to 0. 53685 Vrj if Pg/ P^ 0. 53685 and to Pg if Pg/P^ 

i 


^n4 ^v, n4^n4't/ ^n. 2. 2^2. 2 ^ 9 . > ”^9 9 I 


'p, 2 . 2 " 2.2 -^2 — 

\^out 


2.2) n4^"^s, n4 


(B79) 

7 ) (B80) 

>0. 53685. 
(B81) 

-P7) 

(B82) 

(B83) 
>0. 53685. 

^out) 

(B84) 
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where Pg is equal to 0. 53685 P 2 2 ^out^^2 2 ^out ^out'^^2 2 

>0. 53685, and Pq^j^ is equal to Pg for STRM3=”TRUE*’ and to P,^ otherwise. 


6 " ^13 ['^6» 

(B85) 

^p, 3. 5 " ^13^’^3. 5 ’ 

(B86) 

^p, 7 " ^13 p7» 

(B87) 

'p, 2.2 " ^13^'^2. 2’®^ 

(B88) 


for TURBJ1="TRUE" or TURBJ2=”TRUE" 

Fj ^3 for MIX="TRUE’,STRM3="FALSE” 

Fn 3 + Fj ^4 for MIX="TRUE',STRM3="TRUE" ^ggg^ 

Fni + Fj ^2 for HBPR=’'TRUE",MIX="FALSE",STRM3="FALSE" 

Fni + F ^4 for HPBR="FALSE”,MIX="FALSE",STRM3="FALSE" 

Fni + Fj^ 2 + ^n4 for HBPR-”TRUE", MIX="FALSE",STRM3="TRUE'' 



P 2 = Constant 

(B90) 


T 2 = Constant 

(B91) 


^ ^A^2. l'^2. 1 

^2 1 

< 1 . X .y 

'' 2.1 

(B92) 

(1 

•t 

, (*,1 - »e - '^blc - *bls - *v2><*‘ + ^2 , , 

(B93) 
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= 0 for HBPR=’’TRUE", STRM3="FALSE" 


(B94) 


w 


v 2 


■4 


^2. 1 " ^2.2 




v 2 


/ ^ 2.1 

^ n L 


1 - ^2. -^^blc -^^bls -^% 2 > 

S, 2. 1 


^2. l(^fl - - ^blc - ^bls - ^v2) 


dt +T 


2 . 1 ,i 


Po = 


RAW3T3 


3 V<; 


W 3 


Jq ^ovb - ^bll - '^bl2^‘*‘ -"^3 . 


Jo W„ L 


+ ^ovb+^bll ^^bl2> 


- T, 


|< *c ‘ ■ *ovb - *bll - *bl 2 > j 


dt +T 


3,i 


^3 "^4^'^3, 0) 


^4 = 




W. 




(B95) 

(B96) 

(B97) 

(B98) 

(B99) 

(BlOO) 


(BlOl) 


T,=r_L p 

1 *4 L 




^4 


-T4(Wb+WF-Wti)j 


dt.T4 . 


^5 = 


^A^gTs 


(B102) 


(BIOS) 


'^5 = ^ (*tl +*bll-*t2)*+W5 1 


(B104) 


' Jo -sL ^ 


bllhg - 


- TgCwjj + Wbii - 


‘^*■^'^5,1 (BIOS) 


p 

V, 


(B106) 


6 = ^ ('^t2+'^bl2 -\l)dt+Wg j 


(B107) 


’• ik 


^2^^6+^bl2^3-%1^6 

S,6 


- ^6^2 +^^^2 - 


dt + Tr . 

6,1 


(BIOS) 


p _^A^3,2'^3.2 
3 . 2 " 

^ 3^2 


(B109) 


3.2=^ (^f2 ^'^bls - "^v3^dt +W3 2 ^ i 


(BllO) 
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^f24.2^^bls^2,l-^v3^3.2 

S.3.2 

-'^3.2^^f2+^bls-^v3> 

RAW 3 . 5 T 3.5 


dt +T 


3.2,i 


3.5 


3.5 


W 


>.5 =/n 


t 


3.5 V/o (^3- ^^**""^3.5.1 


(Bill) 

(B112) 

(B113) 


>■5 Jo W 3_5 


*v3‘'3.2 - *n2*'3 .j . - *„2) 


"v,3. 5 


2 . 2 


"^3. 5,i 


^A^ 2 . 2 '^ 2 . 2 
V 2.2 


W 


!.2=/o*<*v2- * 1 . 4 )" + '*'2.2.1 


(B114) 

(B115) 

(B116) 


2.2 


=r — 

Jo W 2_2 


*,21.2 1 - 

^v, 2. 2 




+ '^2.2,1 


^7 = 


RaW^t, 


(B117) 

(B118) 




*7 =i 

't 

(w. + 

0 

*„i - 

*n3>* ■* 

in^in ^ 

^ni^e - 

*n3'>7 

- T,(*i„ 


'nl - "'n3^ 


dt + T 


'v,7 


(B119) 

7,i 

(B120) 
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where is equal to ^^^ 2^3 5 HBPR="TRUE", to w^4^^2. 2 HBPR 

and MK='’TRUE," and to 0 otherwise. 


P 


8 


= Constant 


w „ = Output of fuel control 
r 


^hJo 


(hi - " ^ti, i 




— (h2 ■ ^fl “ ^ i 

^It2 ^ 


Nfi = KfiN^2 


^f2 ■ *^f2^t2 


Nc = Nti 


^v3 " 


/ 


h.2 - ^ 3.5 




v3 


For SUPER ='TRUE" 


prfl = a + b(flnp) 


= Constant = 


"FALSE" 

(B121) 

(B122) 

(B123) 

(B124) 

(B125) 

(B126) 

(B127) 

(B128) 

(B129) 

(B130) 
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(B131) 


+ w 


blc + '^bls 


(B132) 


APPENDIX C 


SCALED- VARIABLE SYMBOLS 


ABLC 

control bleed nozzle area 

ABLS 

inter stream bleed flow coefficient 

AI(i) 

DAC initial condition array, i = 1 to 24 

AN^ 

nozzle area, I = 1 to 4 

C(i) 

analog coefficient array, i = 1 to 57 

CPAB 

specific heat at constant pressure of combustor air 

CP^ 

specific heat at constant pressure at station I 

CVAB 

specific heat at constant volume of combustor air 

CVi 

specific heat at constant volume at station I 

DTQTl 

high-pressure-turbine differential torque 

DTQT2 

low-pressure-turbine differential torque 

DW^ 

stored mass derivative at station I 

FARI 

fuel-air ratio at station I 

FG3 

compressor discharge bleed flow coefficient 

FPC 

high-pressure- compressor map flow parameter 

FPFl 

low-pressure- compressor map flow parameter 

FPF2 

fan map flow parameter 

FIBL 

control bleed flow parameter 

F2BL 

control bleed thrust parameter 

FINI 

nozzle flow parameter, I = 1 to 4 

F2NI 

nozzle thrust parameter, I = 1 to 4 

GMAB 

specific-heat ratio of combustor air 

GMX 

specific-heat ratio at station I 

HAB 

enthalpy of combustor air 

HABS 

rescaled enthalpy of combustor air 



HBPR 

m 

HIP 

HB 

1HR.H 

ISF 

K(k) 

JY(k) 

kx 

LERR 

MIX 

MN 

MODE 

NCV 

NMP 

NPT 

NSC 

NTBL 

NX(k) 

NY(k) 

N1 

N2 

PCNC 

PCNFl 

PCNF2 


logical variable denoting separate characteristics for fan and 
low-pressure compressor (fan hub) 

enthalpy at station I 

enthalpy at inlet to station I volume 

rescaled enthalpy at station I 

hybrid linkage routine error flag 

overflow error flag (tested by LERR) 

initial pressure-ratio search index array for MAPFUN routine, 
k = 1 to 12 

initial speed search index array for MAPFUN routine, k = 1 to 12 
integer specifying location of component data to be scaled 
logical function error detection routine 

logical variable denoting mixing of streams at turbofan discharge 
component map index, MN = 1 to 10 
integer indicating analog mode 
number of curves per map, NCV=8 

number of component maps, NMP< 10 (MAPFUN cai^ble of 
handling NMP < 12) 

number of points per curve, NPT=8 

number of digital coefficients, NSC < 125 

number of map data points stored in VALS array, 
NTBL=NCV*(2*NPT +1) 

number of points per curve array, k = 1 to NMP 

number of curves per map array, K = 1 to NMP 

index on lowest component map 

index on highest component map 

fraction of design high- pressure- compressor speed 

fraction of design low-pressure-compressor speed 

fraction of design fan speed 
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PCNTl 

PCNT2 

PI 

PRC 

PRFl 

PRF2 

PSBL 

PSNT 

SC(i) 

SENSW 

SF^, 

SFx 

SPLIT 

STRM3 

SUPER 

TAB 

TAUNI 

TAUT 

T^ 

TIDN 

T^S 

TRFl 

TRQC 

TRQFT 

TRQFl 

TRQF2 

TRQTl 

TRQT2 


fraction of design high-pressure-turbine speed 
fraction of design low-pressure-turbine speed 
total pressure at station I 
high-pressure- compressor map pressure ratio 
low-pressure-compressor map pressure ratio 
fan map pressure ratio 

static pressure at control bleed nozzle throat 

static pressure at nozzle throat, I = 1 to 4 

digital coefficient array, i = 1 to NSC 

logical function routine for testing sense- switch positions 

scale factor on time t, t' = SF^.t 

scale factor on variable x, X = x/SF^ (appropriate units) 

logical variable denoting fan maps based on total engine flow 

logical variable denoting three-stream turbofan 

logical variable denoting linear characteristic for low-pressure 
compressor 

combustor air temperature 
nozzle thrust, I = 1 to 4 
total thrust 

total temperature at station I 

specific temperature derivative at station I 

rescaled temperature at station I 

low-pressure- compressor ideal temperature rise parameter 
high-pressure-compressor torque 
fan torque based on total inlet flow 
low-pressure-compressor torque 
fan torque 

high-pressure-turbine torque 
low-pressure-turbine torque 



TURBJl 

TURBJ2 

t' 

VALS(i) 

WDB 

WDBLC 

WDBLS 

WDBLl 

WDBL2 

WDC 

WDF 

WDFl 

WDF2 

WDFT 

Wi 

WDOT 

WDOUT 

WDOVB 

WDTl 

WDT2 

WDV2 

WDV3 

X 

XDAC 

XSC 

XVALS (i,j,k) 
YSC 

YVALS (j,k) 


logical variable denoting one-spool turbojet configuration 
logical variable denoting two-spool turbojet configuration 
computer time 

array containing unsealed component data, i = 1 to NTBL 

combustor airflow 

control bleed flow 

interstream bleed flow 

high-pressure-turbine cooling bleed flow 

low-pressure-turbine cooling bleed flow 

high-pressure- compressor flow 

fuel flow 

low-pressure- compressor flow 
fan flow 

total engine inlet flow 

stored mass at station I 

nozzle flow, I = 1 to 4 

flow out of control volume 

overboard bleed flow 

high-pressure-turbine flow 

low-pressure-turbine flow 

second- stream duct flow 

third- stream duct flow 

scaled variable X = x/SF^ 

variable transmitted through extra DAC 

map pressure-ratio scale factor 

scaled map pressure-ratio array, i = 1 to NPT; 
j = 1 to NCV; k = 1 to NMP 

map speed scale factor 

scaled map speed parameter array, j = 1 to NCV; k = 1 to NMP 


86 



zsc 

ZVALS(i,j, 


map output scale factor 

scaled map output array, i = 1 to NPT; j = 1 to NCV; k 


1 to NMP 
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APPENDIX D 


DEFINITIONS OF DIGITAL AND ANALOG COEFFICIENTS 


Digital Coefficients 


SC(1) =-^ 
SFt2 

SC(2) 

^%DN4 


SC(3) = 


SF 


P7 


SF 


P22 


SC(4) = Kfi 


^t2, ctes 
^fl, des 


SC(5) = Kj2 


^t2, des 
^f2, des 


SC(6) = 


^bll^^P3 


Sr 


^®^WDBL1 V ^^T3 ^A 


SC(7) = 


^bl2 ®^P3 


2SF 


WDBL2 


SC(8) = 


2SF, 


®^T3 *^A 


/ 

B Y®^T3 


SC(9) 


SFps 

SFpg 

®^P21 

®^P22 


SC(IO) =!!z2Si£_ 
®^ABLC®^P21 


^l/2qF 4?^/^ <?1? 

SC(ll) - b ^^WDB _ ^%DB 





SC»2) . _1!wDF 

®^WDB®^FAR4 


SC(13) = 


®^WDBL1 


WDTl 


SC(14). =. 


WDBL2 

VdT2 


SC(15) 


SC(16) = 


. ^WDNl 
®^AN1®^P6 


fSF^g 


SC(17) =. 


SC(18) 

®^AN2^^P35 


®^T35 ^A 
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SC(19) 


SC(22) 


^^TAUN2 

®^WDN2^v, n2 T ®^T35^^CP35 

SC(20) =.!!tAUN2_ 

®^P35^^AN2 


SC(21) = 


^^WDN4 

^^AN4®^P22 



^^TAUN4 
^^WDN4^v, n4 




®^T22®^CP22 


SC(23) = 


^^TAUN4 

SFp22^^AN4 


SC(24) = ^‘'tAUNI 
^^TAUT 

SC(25) =!!tATO2 
^^TAUT 

SC(26) =!!l^ 
®^TAUT 

SC(27) 

®^WDN2 


SC(28) = 


SF 

SF 


P7 


SC(29) =. 


SE 


WDN3 


®^AN3®^P7 




SFt7 


SC(30) = 


SF, 


TAUN3 






®^WDN3‘^v, n3 If ®^T7®^CP7 


SC(31) =. 


SF 


TAUN3 


SFp7SFAN3 


SC(32) =. 


SF 


" bis 

ABLS 


SC(33)>!!wM 

®%DN3 

SC(34)=!!«™! 

^%DN3 

SC(35) 

®^WDN3 


SC(36) = 


^^TAUN3 

®^TAUT 


SC(37) =. 


SF, 


TAUNl 




®^WDNl^v, nl V 2J SF,pgSF^pg 


SC(38) = 


SF, 


TAUNl 


SFpgSFAfji 



SC(39) = 


WDBLl 

®^WDC 


WDBL2 

^^WDC 


WDOVB 

^%DC 


®^WDB^%AB 


WDF 

%DB 


^^WDBL1^^T3 

^^WDT1®%5 


. ^^WDBL2^^T3 
®^WDT2^^T6 


SC(47) =. 


WDC 


^^WDFl 


SC(48) = 


WDV2 

%DF1 


SC(49) = 


WDBLC 


^^WDFl 
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SC(57) 


SC(58) 


sc(50) 

®^TRQT2 


SC(51) = ^^trqfz < V 
®^TRQT2 


SC(52) = 


^^TAUN2 

®^TAUT 


SC(53) = 


SF 

SF 


P5 

P4 


SC(54) =. 


SF 


P6 


SF 


P5 


SC(55) = 


^tl.d es^NHPCF 

®^T1NP y^^T4 


SC(56) = 


Nt2,desCNLPCF 

®^T2NP J^^T5 


SFppTiSFp4Nti ®^PCNT1 

SFt5SF^DT2'^P^P^^ 

SFFPT2^^P5^t2,des^^^^^^ ^^PCNT2 


93 



SC(59) =. 


WDFl 


WDFT 


SC(60) = — WDF2 
^^WDFT 


SPLIT = "TRUE' 


SC(61) =. 


SFHPTlNtl,des v/SF^ 4 DHHPCF CNHPCF SFp^^^^ 


lOOOSF 


DHTl 


SC(62) 


S%PT2Nt2,des '/SF^5 DHLPCF CNLPCF SFp^j^2 


lOOOSF 


DHT2 


SC(63) =• 


DHTl 


SC(64) = 

SFt5 


SC(65) = 


1 5 J SFp^,y iSF^p^ ^ 

^®^TRQTl^tl,des®^PCNTl 


SC(66) =. 


15 J SFpjjpgSF^P^g 
^ ®^TRQT2^t2, des^^PCNT2 


SFp2i S^P21 


SC(67) = £-£i ; flfi— for SUPER="TRUE" 

SFpRFiSFp2PRICF 15SFp2 
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SC(68) = 





SC(69) =. 


SF 


P32 


® ^PRF^ Fp2PRFC F 


SC(70) =. 


SF, 


T2 


2T 


2,des 


SC(71) =. 


SF 


T21 


2T 


2, 1, des 


SC(72) =. 


SF, 


T2 


2T 


2, des 


SC(73) = 


'std 


SFfpfiSFp2 


^^WDFl V ”^2, des 


SC(74) = 


^^FPC^^P21 


std 


2. l,des 


sc(^,=1!zzi£!e^ 

S%dF 2* V '^2,deE 

SC (76) = 

SC{77) - /3^ 


*Use SF 


WDFT 


for SPLIT=*'TRUE. " 


WAICF 


WACCF 


WAFCF 
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SC(78) = Pj2 


SC(79) = 


S^P21 

15SFp2 


SC(80) = 


SF 

15SF 


P3 

P21 


^^P32 

SC(81) = 

15SFp2 


SC(82) - ^^TRQFl^fl, des^^PCNFl "" 
®^WDF1®%21S^°'^ 

SC(83) = ^^TRQC^c, des^^PCNC 

g^(34) _ ^^TRQF2*%, des^^PCNF2 ^ 
^^WDF2*®%32S^®*^ 

sc(85) =_^niL_ 

®^AN4 

SC(86) — 

®^ABLC 

SF 

SC(87) = — for SPUT="TRUE" 
®^TRQFT 


*Use SF^ppp and SFppgpp for SPLIT="TRUE. ” 
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SC(88) = 


^ v3^^^WDV3 
|/®^P32 


^v3^^^WDV3 
|/ ®^P35 


SC(89) = 


®^WDV2 


^1/2 op 

^v2 ^*WDV2 


/ 


SF 


P21 


/ 


SF 


P22 


SC(90) - r) 


sc 


gC(9i) ^ ^^TRQFl^fl. des^^PCNFl 
30J SFjj2iSFwj)fi 


for SUPER=”TRUE” 


1 - 


SC(92) = 


PRICE 


^^PRFl 


1 


SC(93) = 


SF 


PRCCF 

PRC 


SC(94) = mFCF 

®^PRF2 


SC(95) 

ETAICF 


SC(96) 

ETACCF 
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SC(97) =_2i25 — 
ETAFCF 


SC(98) =. 


SF 


P32 


SF 


P21 


SC(99) = 


SF, 


WDBLS 


^ SFy-,<j^SF 


P21^'ABLS 


SP(inn) ^ ^^WDBLS _ ^^WDBLS 

®^WDF2 ^^WDV3 


SC(lOl) - ^^WDBI^ _ ^^WDBLS 
^%DF1 ^^WDC 


SC(102) =• 


*^n3 


SF 


AN3 


SC(103) = 


SF 


P8 


Analog Coefficients 


C(l) = 


^%DF1 


98 


Wo 1 i 

C(2) =. 


SF, 


W21 


C(3)=!!wDfl_ 


C(4) = 


_ '^2. 1, i 


SF, 


T21 


C(5) = 


^%DV2 

®^W22®^t' 


Wo 2 i 

C(6) = — ±Jl± 


SF 


W22 


C(7) = ^^W _ Dy 2_ 
®^W22^^t' 


C(8) = 


'^2.2, i 

SFt22 


C(9) = 




SF 


W3 


99 



C(10) = — §Li 

SFt3 


C(ll) = 


SF, 


WDC 


SF^gSFj, (10) 


C(12) = ffwDC 


C(13) = 


W 


4, i 


SF 


W4 


C(14) = 




SF, 


T4 


C(15) = 


^^WDB 
SFw4SFt, (10) 


C(16) = 


^^WDB 
SF^^SFj, (10) 


C(17) = 


W 


5, i 


SF, 


W5 
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T 

C(18) 


SF, 


T5 


C(19) =_!!w5T1_ 

sf'wssr'i. (10) 


C(20) = .... ! .!wDTj_ 

SF^jSFj, (10) 


C(21) 





C(22) = 


^6 

SF 


T6 


C(23) = 


^%DN1 
SF^gSF^, (10) 


C(24) = 


^%DN1 
SF^gSF^, (10) 


C(25) =!!w2X5_ 
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C(26) = 


^A^^W21^%21 

SFp2lV2. 1 (10)' 


5 i 

C(27) = — 


SF, 


W35 


C(28) = 


^A^^W35®^T35 


^^P35^3. 5 


( 10 ) 


C(29) = 


SF 


WDF2 


^^W32®^t' 


C(30) = 


^A®^W22^^T22 


®^P22^2. 2 


( 10 ) 


C(31) = 


W 


3. 2, i 


SF 


W32 


cr 32 ) = ^A^^W32^^T32 
' SFP32V3.2C0) 


C(33) = 


SF 


WDN3 




C(34) = ^A^^W3^^T3 
SFpgVg (10) 
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C(35) = 


W 


7 , i 


SF, 


W7 


C(36) = _ ^A^^W7^^T7 
SFp^V^ ( 10 ) 


C(37) = 


30SFtrqti 

des^^t'^^PCNTl 


C(38) = 

SFP 4 V 4 ( 10 ) 


C(39) = 


N, 


tl, i 


Ntl, desSF 


PCNTl 


C(40) = 


30SFtrqt2 

’^^12^12, des^^t'^^PCNT2 


C(41) 

SFpgVg ( 10 ) 


C(42) = 


N 


12, i 


^12, des®^PCNT2 
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C(44) = 


C(45) 


C(46) = 


C(47) 


C(48) = 


^^WDV3 


'^3. 5, i 
®%35 


^%DF2 

S^W32®^t’ 

_ '^3. 2, i 

SFt32 


®%DN3 



C(53) = 


^ TT1 NP 
C(50) = 

®^PRF1 


C(51) =, ^ 

®^PRF1 


C(52) 

(10)SFp2i 




SUPER=’'TRUE" 
I a + b(flnp) = prfl 


^fl^t2, des 
^fl, des 



10 


C(54) =. 


n2 


SF 


AN2 


C(55) = 



C(56) =. 


Si 


SF 


ANl 


w. 


C(57) =. 


Ei± 


SF 


WDF 
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APPENDIX E 


COMPUTER LISTINGS OF DATA INPUT PROGRAM, 


MAIN PROGRAM, AND SUBROUTINES 


0EMERALIZED TURBOFAN ENGINE DATA INPUT PROOTAM 


DIMENSION VALS(ISS) 

SCALED FRACTION SC ( 125) ,AI (24) ,XVALS(I 1,8, I2>,YVALS(8, 12) , 

1 ZVALS(1§,8,12) 

CWM0N/MAPS/XVALS,YVALS,ZVALS,IX(12),JY(I2),NX(12),NY(12) 

C OMMON/ AZ /SC /BZ /MN/CZ / AI 
DEFINE ISFC635) 

LOGICAL SENSW,LERR 

C*****DETERMI NE MAP NUMBERS AS USED IN MAIN PROGRAM 

1 TYPE 2 

2 FORMAT (/3X,S#HDEPRESS SENSE SWITCHES TO SELECT CONFIGURATION. THE 
IN R-S-R./) 

PAUSE 

5 IF(SENSW(5)) GO TO 7 
IF(SENSW(2)) GO TO 8 
N1 = I 
N2=8 

GO TO II 

7 NI = 3 
N2r« 

GO TO 11 

8 NZ = 1 1 

1F(SENSW(7)) GO TO 9 
Nl = l 

60 TO If 

9 Nl=3 

II NMP=N2>NK1 
ISF=I 
TYPE 13 

13 FORMAT (/3X,S3HPLACE DATA TAPE FOR DIGITAL COEFFICIENTS IN HSPTR. 

I THEN R-S-R,/) 

PAUSE 1 

READ (4,14)NSC 

14 F0RMAT(I4) 

READ (4,I9)(SC(I),I=1 ,NSC) 

19 F0RMAT(9S8) 

TYPE 16 

16 F0RMAT(/3X,65HPLACE DATA TAPE FOR DAC INITIAL CONDITIONS IN HSPTR. 
1 THEN R-S-R./) 

PAUSE 2 

READ <4,15)(AI(I),I=1 ,24) 

C*<n»t<niiimTIALlZE FOR READING MAP DATA 
NCVr8 
NPT=8 

NTBL=NCV*(2*NPT+1) 
type |7,N1,N2 

17 FORMAT (/3X, 3 IHPLACE DATA TAPES FOR MAPS N0. ,12, 1H-,I2, 1 IH IN HSPT 
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IR*/) 

TYPE 18 

18 F0RMAT(/3X,42HDEPRESS SSW(A) FOR DATA LIST* THER R-S-R,/> 
DO 28 NrNI,N2 

PAUSE 3 

READ (4,I9)XSC,YSC,ZSC 

19 F0RMAT(3F9.2) 

C***4>4fiEAD Y VALUES 

READ (4,2IHVALS(I),I = 1 ,NCV) 

2# FflRMAT<8F9,3> 

J=NCVfI 

C**4ni»l>READ X,Z PAIRS 

READ <4,21)(VALS(I),I=J,NTBL) 

21 F0RMAT(8F9,4) 

C<t***»scALE HAP VALUES 

DO 22 Jrl,NCV 
nn 99 t MPT 

IF<VALS’(J)*EQ.YSC) QO TO 31 
YVALS ( J , W )r VALS ( J ) /YSC 
QO TO 31 

31 YVALS (J*R>r.99999S 

31 KXsMCVM-24«(<J-l)*NPT+I)-l 
IF(VALS(KX>.EQ,XSC> 60 TO 32 
XVALS(I,J,N)rVALS(KX)/XSC 

60 TO 33 

32 XVALSn,J,M)r,99995S 

33 IF<VALS(KX+n,EQ,ZSC) 60 TO 34 
ZVALS(I,J,N)rVALS(KX+l)/ZSC 

60 TO 22 

34 ZVALS (I *J*N)r ,99999s 

22 CONTINUE 

C4c4nMt*TEST FOR SCALED FRACTION OVERFLOW 
1F(.N0T.LERR(I3)) 60 TO 24 
TYPE 23 

23 F0RMAT(24HSCALED FRACTION OVERFLOW) 

PAUSE 4 

ISFsl 

24 lX<N)r| 

JY<N)=I 
NX(H)=MPT 
NY(N)s NCV 

IF(.HOT*SENSW(I)> 60 TO 26 
TYPE 25,(YVALS(J,N),Jrl ,NCV) 

29 FCRHAT(8S7) 

DO 26 Jr|,NCV 

TYPE 25 , (XVALS < I , J , N) ,Z VALS ( I , J , N) , 1= I , NPT) 

26 CONTINUE 
TYPE 27 

27 FCRHAT(/3X,I9HDATA IS LOADED,/) 

PAUSE 9 

END 
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GENERALIZED TURBOJET OR TURBOFAR ENGINE SIHULATION 


C««*4»itADC VARIABLES . - 

SCALED FRACTION P2I ,T21 ,P22,T22,P5,T3,P4,T4,P5,T9,PG,T6, 

I P59,T39,P32.T32,P7,T7,PCNT1,PCNT2,AN2,P2,ANI,WDF 

C^4t#4E4cDAC VARIABLES • 

SCALED FRACTION T2IDN,DW21 ,T220N,DW22,T390N,DW39,T32DN,DV32, 

1 T3DN,DV3,T40N,DW4,T5DN,DW9,T60N,DW€,T7DN,DW7,0T0T1,DTQT2, 

2 PRC.FPC.TAUT.XDAC 

SCALED FRACTION ABLS,T2,P8,SC(I29) ,PCNFI ,PCNF2,PCNC,T2IS, 

1 T3S,T9S,T6S,T35S,T7S,CP21 ,QM2I ,CV2I ,CP39,CV35,6M39,CP3,CV3, 

2 GM3,FG3,WDBL1 ,VDBL2,WD0VB,PSN4,F2N4,FI N4,WDBLC,VDB, 

3 FAR4»CP4,CV4,QM4,WDFI ,TRQF1 ,VDC,TRQC,VDTI ,TRQTI ,FAR9, 

4 CP5,CV5,0M5,WDT2,TRQT2,FARS,CP6,CV«,QW6,PSNI ,F2N1,FINI , 

5 WDN1,VDF2,H3,H3S,H2I ,H2!S,H4,H4S,H21P,H3P,M5P,H6P,H9, 

6 TRQP2,PSN2,F2N2,FIN2,WDN2,WDM4,FAR7,CP7,CV7,QM7,WDF^TRQFT, 

7 PSN3,F2N3,F1N3,WDN3,TAB,CPAB,CWAB,CVAB,PRF1,FPFI,PRF2,FPF2, 
g PRFIS,AI(24>,XVALS(1 I,8,12),YVALS(8,12>,ZVALS(I«,8*I2),SSQRT 

SCALED FRACTION PSBL,F2BL,F!BL,AN4, ABLC,TAUNI ,TAUN4, 

1 TAUN2.TAUN3,T52S,T22S,CP22,CV22,GM22,CP32,CV32,GM32,WDV2, 

2 ¥DV3,WbOUT,TRFl ,H5S,HG,H8S,H32,H32S,H39 ,H39S,H7,H7S, 

3 H22,H22S.HAB,HABS,T2S,CP2,CV2,QM2,H2S,M2,H32P,AN3,VDBLS 
C OMMON/MAPS/XVALS , YVALS , Z VALS, IX(I2),JY(12),NX<I2),NY(I2) 

C OMMON/AZ /SC /BZ /MN/CZ / AI 

DEFINE ISF( *635) 

LOGICAL SENSW,LERR,HBPR,MIX,STRM3,TURBJI,T0RBJ2, SUPER, SPLIT 
C*««4t4riNITIALIZE ANALOG 

CALL QSHYlN(IERR,68f) 

CALL QSRUNdERR) 

CALL QSSECN(IERR) 

C<t*>w*ESTABLISH CONFIGURATION 
TYPE 17 

17 F0RWAT(/3X,6fHDEPRESS SENSE SNITCHES TO SELECT CONFIGURATION. 
IN R-S«R,/> 

PAUSE 

IF(SENSN<2)) GO TO I 
HBPRs, FALSE, 

CALL QNCLL(2,.FALSE.,IERR) 

GO TO 2 


1 HBPRs.TRUE. 

CALL QVCLL(2,.TRUE.,IERR) 

2 IF(SENSW(3)) GO TO 3 
«1X= .FALSE, 

CALL QNCLL(3,.FALSE.,IERR) 
GO TO 4 


3 MiXs.TRUE. 

CALL QWCLL<3,,TRUE.,IERR) 

4 IF(SENSW<4>) GO TO 9 
STR«3=. FALSE, 

CALL QWCLL(4..FALSE,,IERR) 
GO TO S 

9 STRHSe.TRUE. 

CALL QVCLL(4,.TRUE.,IERR) 

5 IF(SENSW(9>) GO TO 7 
TURBJlr. FALSE. 

CALL 0WCLL<9,,FALSE,,IERR) 


THE 
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GO TO g 

7 TURBJ1=.TRUE. 

CALL QWCLL(5,.TRUE.,IERR) 
g 1F<SENSW(€>) GO TO 9 
TURBJZr. FALSE. 

CALL QWCLLCS,. FALSE., lERR) 

60 TO 10 
9 TURBJ2=.TRUE. 

CALL QWCL^(6,.TRUE.,IERR) 

If IF(SENSW(7)) GO TO II 
SUPERr, FALSE. 

CALL QWCLL(7,. FALSE., lERR) 

60 TO 12 

11 SUPERS. TRUE. 

CALL QWCLL(7,.TRUE.,IERR) 

12 IF(SENSW(g)) 60 TO 13 
SPLITS. FALSE. 

60 TO 14 
15 SPLITS TRUE 

C»****gTOP PROGRAM IF COMFIGURATIOW SPECIFIED IS MOT ALLOWED 
14 IF(TURBJI .AND.TURBJ2) STOPl 

IF((TURBJI .0R.TURBJ2).AND.(HBPR.0R.MIX.0R.STRMS.0R.SUPER.0R. 
1 SPLIT)) ST0P2 

I F( . ROT . HBPR . AND . (STRMS . OR .SUPER . OR .SPLIT) ) STOPS 
IF<STRM3.AND.SUPER) ST0P4 

c*****initialize required input /output 

PgsSCdfS) 

T2sSC(l) 

TAUNls.flff IS 

TAUN2s.f|fffS 

TAUNSs.llfllS 

TAUN4s.ffffllS 

WDFls.ffffffS 

WDFgs.llflfg 

WDNls.lflffg 

WDN2s.fllflS 

WDNSs.IMfIS 

VDN4s.ff||fg 

AN4sSC(g9) 

AN3sSC(lf2) 

ABLCsSC(gg) 

ABLSsSC(S2) 

WDV2:.lllffS 

WDVSs.lfffIfS 

TRQFgs.ffllfS 

TRQFls.lMIfg 

TRQT2s.ffflfflS 

WDBLCs.fffffS 

WDBLSs.ffflllS 

19 CALL QWBDAS(AI,f ,24,IERR) 

21 CALL QSTDA 
31 CALL QSIC(IERR) 

ISF=f 
TYPE 39 

39 F0RMAT(/5X,44HPR0CEED TO DYNAMIC PART OF PROGRAM BY R-S-R./) 
PAUSE 

4f CALL QRBADS(P21 ,I,24,IERR) 

42 PCNF|sSC(4)*PCNT2 



PCNF2=SC(5)*PCNT2 

PCMCrPCNTI 

T2S= ,2ilMS*T2 

T22Sr.2HilS*T22 

T52S=,2f#HS*T52 

T35S=,2*lifS*T35 

T3Ss.4lilfS*T3 

T5S=.8iiMS*T5 

T6S=.5fiilS*T6 

T7Ss*5llliS*T7 

CALL PR0C0M(T2S,,ilfliS,CP2,CV2,QM2,H2S) 

H2sH2S/.33333S 

CALL PR0C0M(T3S, .iilliS,CP3,CV3,QM3,H3S) 

H3=H3S/.66«67S 

FQ3s.A42fl8S+,34S74S*QM3 ...... 

WDBLI = (<SC(S)*FQ3*P3)/SSQRT(T3))/.5MMS 
VD6L2= ( (SC (7)*FQ3*P3)/SSQRT(T3) )/,5IIMS 
WDOVB= ((SC (8)*FQ3*P3)/SSQRT(T3))/«9fAliS 

49 IF (.ROT, SUPER) 60 TO 99 
PRFISr (SC (67)*P2 1 ) /P2 
CALL TRAT(PRFlS,.7liffS,TRFI) 

T2l = ((»29MiS*TRFl/SC(9l)+.2it4*S)<'T2)/.2flffS 

CALL"pR0C0M(T2?S,,fMiiS,CP2!,CV2! ,6M2I,H2IS) 

H2! = MeiS/.33333S 
IF(TURBJI> 60 TO 9i 

CALL «0ZZL(P21,P8,SC(9),F1BLjF2BL,MBL) 

WDBLCr((P2 l*FIBL*ABLC)/SSQRT(T2n>/SC(in 
91 WDBs(SSQRT(PS-P4))/SC(11) 

FAR4= (SC ( 1 2)*W0F) /WDB 

CALL PR0C0«(T4,FAR4,CP4,CV4,8M4,H4S) 

H4=H4S*.6liflS 

MR=I 

IF(TURBJI .OR.SUPER) MR=3 

CALL**FRCMP(1 ,P2,T2,H2,P2I ,T2I ,PCHF1 ,WDF1 ,H2IP,TRQFI ,PRF1 ,FPF1) 

99 CALl"fRCMP(2,P2I,T21.H21,P3,T3,PCNC,WDC,H3P,TRQC,PRC,FPC) 

IF(,MOT.HBPR) 60 TO 96 

lF(ABLS,LE,,liiilS) 60 TO 96 

iwLSr ( ABLS*SSQRT(P2 1 -SC (98)*P32) )/SC (99) 

96 1F(. ROT .SUPER) 60 TO 97 
iDF!"SDC*SC(47)+WDBLC*SC(49)+WDBLS*SC(llI) 

TRQF1=(((H2 I-H2)/SC(91))*WDF1 )/PCRFI 

97 riR=MR+l 

CALl'tURB(!.P4.T4,H4.P9.PCNTI,WDTI.H9P,TRQTI) 

FAR9=(FAR4*.9lfMS)/(.9lff«S+(.9«fi#S+.f29liS*FAR4)*SC(!3)* 

*CAU**PR0C0M(T9S,FAR9,CP9»CV9,6M9,H9S) 

M9=H9S*.79#MS 
IF(.ROT.TURBJl) 60 TO !•! 

T6=T9/.629flS 

T6S=,9MliS*T6 

FAR6=FAR9 

P6=P9 

GO TO III 
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Iff CALL TURB(2,P5,T5,H5,P6,PCNT2,WDT2,H6P.TRQT2) 

FAR6J(FAR5^.5ffflS)/(.5lflffS+<.5ffffS+.f25ffS*FAR9)*SC(I4>* 

I WDBL2/WDT2) 

Ilf CALL PR0C0M(T6S,FAR6,CP6,CV6,GM6,H6S) 

H6=M6S/.B3533S 

CALL IIOZZL<P6,P7,SC<15),FlNl,F2Nl,PSHn 
WBRls({P6*FI Nl*AWl)/SSQRT(T6))/SC(!€) 

CALL**PR0C0m”t 32S JfffffS,CP32,CV32,GM32,H32S> 

CALL”pRW0M<T35l,.ff#ffS,CP35,CV39,QM39,H39S) 

CALL^FIiCMP<3,P2 tT2«H2 ,P32,T32,PCNF2,WDF2tH32P,TRQF2,PRF2,FPF2) 

IF(,HOT. SPLIT) GO TO 19f 

WDFT=WF2 

TRQFT=TR0F2 

WDF2= ( WBFT-SC ( 99 )*WOFI ) /SC ( 6 f ) 
TRQF2s(SC(6f)<«TRQFT*WDF2/WDFT)/SC<87) 

S^^S(P3SP^sc(l7^^^!^,^2»^.Ps;^^ 

WII2= ( (P39*F 1 N2*AM2> /SSQRT (T 39 ) ) /SC < 1 8) 

19f WDV3=<SSQRT(P32-P39))/SC(88) 

DW39=WDV5-WDN2 

T39DR=(H32*WDV3-M39*WOM2)/CV39-T39<«'DV39 

TAURi=(W0Nl*F2Hl«SQRT(CP6*T6))/SC(57)+(AMl*CPSRI-SC(l9)* 

*TAU^=(WDM2iF2R2*SSQRT<CP39*T39))/SC(19)+(All2*(PSM2-SC(l7)* 

I P7»/SC<2f) 

IF(STRM3) GO TO 21 f 

TAUT=SC (24)<'TAUIil+SC<92)*TAUR2 

2ff FAR7r(FAR«*,9ffffS)/(.9ffffS+<.9ffffS+.f29ffS*FAR6)/(SC(27)* 

I WDRI/WDN2)) 

CALL PROCOW (T7S,FAR7,CP7,CV7,QM7,H7S) 

H7rH7S/*83333S 

CALL »OZZL(P7,P8,SC(28)fFI»l3,F2tl3,PSN3) 
WDN3s<(P7*F|R3*AH3)/SSQRT(T7))/SC<29) 

TAUR3=(W0N3<tF2ll3*SSQRT(CP7*T7))/SC(3f)+<AN3*(PSM3-SC<28)>*' 

I P8))/SC(3I) 

2f6 DW7s-¥DN3+SC(34)>«rt>NI+SC(39)*WDII2 

T7DH2 ( *4f f f fS*SC (39)*H39*V#DN2+SC (34)*M8*WDRI-H74iWDN3) 

1 /CW7-T7*DW7 
1F(STRM3> GO TO 21 f 
TAWT=TAUN3 
60 TO 228 

21 f CALL H0ZZL(P22,P8,SC(9),FI(I4,F2N4,PSII4) 

211 CALL PR0C0W(T22s/.ffflfS.CP22.CV22.6M22,H2BS) 

H22=1C2S/.53333S 

WDII4r ( (P22*FI N4*AN4 > /SSORT (T22 > ) /SC <21 ) 

WO V2= (SSQR T ( P2 1 - P2 2 ) ) /SC < 89 ) 

DW22:WDV2-WDN4 

T22DN= <H2 I*WDV2-H22*W0N4) /CV22-T22*0^2 

IF( NOT HBPR) 60 TO 219 

TAUN4=(WDN4*F2M4*SSQRT(CP22*T22))/SC<22)+<AM4*<PSH4-SC(9)* 
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! P8))/SC(23) 

IF(MIX) 00 TO 212 

TAUT=SC(24>*TAUM1+SC<25)*TAII»2+SC(26)*TAUM4 


GO TO 227 

212 TAUT=SC (5«)*TAUN3+SC (26)*TAUN4 
00 TO 227 

215 1F(TURBJ1.0R.TURBJ2) GO TO 22A 

CALL N0ZZL(P22,P7,SC(3),FIM4,F2N4,PSN4) 

GO TO 21 1 

219 IF(MIX> GO TO 222 ^ , 

TAUN4= (WDIi4*F2N4*SSQRT <CP22*T22) ) /SC (22)+( AN4*(PSM4-SC (3)* 


I P7))/SC(23) 

221 TAUI»1:(WDN1*F2M1*SSQRT(CP6*T6)7/SC(37)+<AN1*(PSN1-SC(15>* 
1 P7))/SC(38> 

TA«T=SC<24>*TAU1I1+SC(26)*TAUN4 


222 FAR7=(FAR6*.5ili#S)/(.9MMS+<«5ffilS+.l25ffS*FAR6)/(SC(2)* 

1 WDK1/VDN4)) 

CALL PROCOM(T7S,FAR7,CP7,CV7,0W7,H7S) 

tl § fl X ^ ^ 

CALL »0ZZL(P7,P8,SC(28),F1N3,F2N3,PS»3) 

WDfl3= ( <P7*F1 N3*AN3) /SSQRT(T7> )/SC <29) 

TAUWr n#DN3*F2ll3*SSQRT(CP7*T7) ) /SC (3l)+< AM3*(PSW3-SC<28)* 

I P8))/SC(31) 

DW7=-W0R3+SC (34)*W0N1 +SC (33 )*WDN4 

T7DM= (.4111 f S*SC <33 )*H22*WDM4+SC <34>*HS*WDm -H7*WDN3) /CV7-T7*DW7 
TAUT=TAUH3 

227 IF(TURBjn 90 TO 23i 
IF(TURBJ2) 60 TO 228 
WD0UT=SC(48>*W0V2 
GO TO 229 

228 VDOUT=.ffMMS 

229 IF(SOPER) GO TO 25# 

0 W 21 = WDFl-SC(47)*WDC-WD00T-SC(49>*WDBLC-SC(lin*WDBLS 

T2 1 DM2 (H21P*WDF1-H21*(SC(47) ♦WDC+WDOUT+SC ( 49 )*WDBLC+SC < 1 # 1 ) ♦ 

I WDBLS>>/CV21-T2i*DW21 

23# DW5rWDC-WDB-SC(59)*W0BLl-SC<4# )*W0BL2-SC(41 )*WDOVB 

TSDMrdlSPit' WDC - H3* ( WDB+S C ( 3 9 ) *WDBL 1 ♦SC ( 4 # ) ♦WDBL2+SC C4 1 ) 

I *WD0VB)>/CV3-T3*DW3 
TABr ( ,4##iiS*T3-T4)*SC (42)+T4 
CALL PROCOM(TAB,.fi#ifS,CPAB,CVAB,QMAB,HABS) 

HAB2HABS4>.Sf«f#S 

DW4rWDB-WDTl+SC(44)*WDF 

T4DM2 (HAB*WDB+SC (43 )#WDF-H4*WDT 1 ) /C V4-T4*DW4 

IF(TURBJl) WDT2rWDMl 

DM5r WOT 1 - W0T2+SC (39)*WDBL 1 

T5DM2 (H5P*VDT l+SC (45>*H3^WDBL 1 - H5*WDT2> /CV5-T5*DW5 

IF(TURBJl) GO TO 235 

DW62WDT2-WDM1-»-SC(4#>*WDBL2 

T6DM2(H6P*WDT2+SC(46)*H3*WDBL2-H6*WDMn/CV6-T«*DW6 
235 OTQTlrTRQTl-TRQC 

DTQT2=TRQT2-SC(5#>*TRQFl-SC(5n*TRGF2 
C^n|,***TFST FOR SCALED FRACTIOM OVERFLOW 
IF(,M0T.LERR(15)) GO TO 251 
C*#>Mr*TEST FOR OPERATE MODE 
CALL QRAMKMODE) 

IF(M0DE,EQ,4) CALL QSH(IERR) 

TYPE 24# 
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24f F0RMAT(24HSCALED FRACTION OVERFLOW) 

PAUSE 

ISFrfl 

IF(MODE.EQ.4) CALL QSOP(IERR) 

25f IF(,NOT.SENSW(l )) GO TO 37f 
C*****TYPE current VALUES 
C*****TEST FOR OPERATE MODE 
CALL OR AMI (MODE) 

IF(M0DE.EQ,4) CALL QSH(IERR) 

261 TYPE 271 

27« F0RMAT(5X,43H WDFl WDC WDB WDTl WDT2 ) 

TYPE 28i, WDFl, WDC, WDB, WDT1,WDT2 

28f F0RMAT(3X,5(S7,2X)/) 

TYPE 29 f 

29« F0RMAT(5X,43H WDNI WDN2 WDN5 WDN4 WDF2 ) 

TYPE 28f,WDNI,WDN2,WDN3,WDN4,WDF2 
TYPE 3ff 

30i F0RMAT(5X,43H P2 P21 P35 P3 P4 ) 

TYPE 28I,P2,P2I,P35,P3,P4 
TYPE 31 » 

311 F0RMAT(5X,43H P5 P6 P7 P22 P32 ) 

TYPE 28f ,P9,P6,P7,P22,P32 
TYPE 32 ff 

320 F0RMAT(5X,43H T35 T3 T4 T5 T6 ) 

TYPE 280,T35,T3,T4,T5,T6 
TYPE 330 

330 P0RMAT(5X,43H T7 T2 T2I TAUNI TAUN2) 

TYPE 280, T7, T2, T21 , TAUNI, TAUN2 
TYPE 340 

340 F0RMAT(5X,43H TAUN3 WDF PCNT! PCNT2 PCNFl) 

TYPE 280,TAUN3,WDF,PCNT1 ,PCNT2, PCNFl 
TYPE 350 

350 F0RMAT(5X,43H PCNF2 PCNC ANl AN2 AN3 ) 

TYPE 280,PCNF2,PCNC,AN1 ,AN2,AN3 
TYPE 360 

360 F0RMAT(5X,43H AN4 ABLC ABLS TAUN4 TAUT ) 

TYPE 280, ANA, ABLC, ABLS, TAUN4, TAUT 
IF(MODE.EQ,4) CALL QSOP(IERR) 

C<n|.**<tOUTPUT TO ANALOG AND RETURN 

370 CALL QVBDAS(T21DN,0,24,IERR) 

CALL CiSfDA r’ 

GO TO 4# - 

END 


GENERALIZED FAN AND COMPRESSOR SUBROUTINE VERSION II 


SUBROUTINE FNCMP(N, PIN, TIN, HI N, POUT, TOUT, SPD, FLO, HOUTP,TORQ,RR,FP) 
SCALED FRACTION PIN, TIN, POUT, TOUT, SPD, FLO, TOUTP,CPOUT, 

1 CVOUT,GMOUT,TORQ,PR,FP,K,TTERM,SSQRT,NP,MAPFUN,TAV,CP, 

2 CV,W,TR,TOUTS,TINS,SC(125),XVALS(10,8,12),YVALS(8,12), 

3 ZVAtS(10,8,12),EFF,HlN,HOUTP,HS,HOUTS,KH,PRS 
COMM6 n/MAPS/XVALS,YVALS,ZVALS,IX( 12),JY(12),NXC12),NY(12) 
COMMON/AZ/SC/BZ/MN 
PRr(SC(N+66)>ttP0UT)/PIN+SC(N+9I) 

TTERM=.70711S*SSQRT(SC(N+69)*TIN) 


I? 
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MPr ( ,?ill#IS*SPD) /TTERM 
FP=MAPFUN(MM,PR,NP) 


MM=MN+I 


EFF=MAPFUN(MN,PR,NP) 

FLO= (SC ( »H-72 )*FP*P I N) /TTERN 
IF(N.EQ,2) GO TO I«f 
K=.2ffllS 


KH=.35553S 
TINS=TIN 
60 TO 210 
!#A K=.4i#flfS 
KHr,<?6667S 
TIMS=.5llf«S*TIN 

299 TAV=K*{T0UT+SC(M+75)*<TINS“T0UT75 

CALL PR0C0M(TAV, .§H««iS,CP,CV,6M,HS) 

PRS= (SC ( M+7g)*P0UT) /PI N 
CALL TRAT(PRS,6M,TR) 

TOUTP= ((SC(N+94)*TR/EFF+.2f»liS)*TIMS)/,2lflAfS 


TOUTS=K*TOUTP 

CALL PR OCOM (TOUTS,, 

TORQ= ( (FLO*(HOUTS- , 

HOUTP=HOUTS/KH 

RETURN 

END 


PliifS,CPOUT,CVOUT,OHOUT,HOUTS> 

33333S*HIN))/SPD)/SC(N+gn 


GENERALIZED TURBINE SUBROUTINE VERSION II 


SUBROUTINE TURB (N, PIN, TIN, HIN, POUT. SPD,^0,H0UTP,T0RQ> 

SCALED FRACTION PIN,TIN,HlN,POUT,SPD,FLO,HOUTP,TORQ,PR,FP,HP,K! 
1 TTERM , SSQRT , NP . MAPFUN , DHT , SC ( 1 25) ,XVALS (11,8,12) , YVALS ( 8,12), 

%OMMON/M^S/XVALS.YVALS.ZVALS.IX(I2),JY(12).NX(12), 

C OMMON/ AZ /SC /BZ /MN 
TTERM=SSQRT(TIN) 

PRr (SC ( N+52 )*POUT ) /PI N 

NP= (SC(N+54)*SPD)/TTERM/,5fll#S 

FP= MAPFUN (MN, PR, NP) 


MN=MN+1 

HP= MAPFUN (MN, PR, NP) 
FL0=((FP*PIN*SPD)/TIN)/SC(N+5S) 
IF(N,NE,1) GO TO 10i 
K|r.8f0IIS 
60 TO 2f# 

I II KI=.625IIS 

211 DHT=(HP*SPD*TTERM)/SC(N+6I) 

HOUTP= (HI N-SC (N+62)*DHT) /K 1 
T0RQ=((SC(N+64)*DHT*FL0)/SPD)/,5IIIIS 


RETURN 

END 


* 
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GENERALIZED CONVERGENT NOZZLE SUBROUTINE VERSION IT 


SUBROUTINE NOZZL (PI N,POUT , A, FI ,F2,PS) 
SCALED FRACTION PIN, POUT , A,F1 ,F2 ,PS,PR 
INTEGER ERRFLG 
PR:(A*POUT)/PIN 
IF(PR,LE,.53685S) GO TO 1« 

PS=A*POUT 
GO TO 2f 
If PRr ,93685s 
PS=PR*PIN 

29 CALL NRMPN (PR, PR N, ERRFLG) 
IF(ERRFLQ.NE.t ) GO TO 3i 
CALL CLCFN(PRN,F2) 

CALL CLCWN(Fl) 

60 TO 5i 
3f TYPE 4f 

4f FORMAT (26HN0ZZLE INPUTS OUT OF RANGE) 
5# RETURN 
END 


PROCOM - CP, CV, GAMMA, « H CALCULATIONS 


SUBROUTINE PROCOM(T,FA,CP,CV,QAM,H) 

SCALED FRACTION T,FA,CP,CV,6AM.H,CPA,HA, AMW,R,TD,CPF,HF 
IF(T,QE,,46flfS) 60 TO 9f 
IF(T,QE,.24lifS) 60 TO 4§ 

CPA= ,48f 68S- ( , 124S4S-T/.98246S)*T 
HAr.fll76S-f(,56558S+,14«79S*T)*T 


GO TO 61 

41 CPA=.4§528S+(,52945S-.38182S*T)*T 
HAr ,»ll 19S+( .56298S+, I 61 5IS*T)*T 
GO TO 61 

51 CPA=,46i63S+(,3»f24S-.15378S*T)*T 
HA=-,fl87IS+(.64968S-»-.l67fflS*T)*T 
61 AMW:,57940S-.II#90S*FA 
R=.f7945S/AMW 
TD=.70MfS-T 

CPF=,93330S-(,29350S+,8I75IS*TD)*TD 

HFr-,f33l5S+(,63624S+,38625S*T)*T 

M= ( • 66667S*HA+ . 1 6667S*HF*FA ) / ( .5 # M 0 S+ . f 25 • 0S*FA ) 

CP= (. 80000S*C PA+, 080 00S>t'CPF*FA)/ (.800005+, 0400 0S»FA) 


CV=CP-R 


6AM=.500 00S>«;P/CV 

RETURN 

END 



C*»***TWO-SPOOL ENGINE ISENTROPIC TEHP RATIO CALCULATIONS 


SUBROUTINE TRATCPRC ,GAM,TR) 

SCALED FRACTION GAM,TR,PRC,S,C,TRC 
INTEGER ERRFL6 
REAL XN 

151 IF(QAM,GE.,G75IIS) GO TO 5il 
IF(PRC,GE,.53333S) GO TO 
S= ,5f526S+,77229S*PRC- ,951 64S*PRC*PRC 
60 TO 711 

2fl S=.58745S+.2385IS*PRC-.i9fl«S*PRC*PRC 

60 TO 7#i 

3fl IF(PRC,GE,,33333S) GO TO 411 

S=,459liS-f,89475S*PRC-(PRC*PRC)/.92895S 
60 TO 7M 

4#f S=,55215S+.2951 9S*PRC- ,1191 9S*PRC*PRC 

7fl C=((,G75IIS-QAM)*S)/,i25flS 
CALL NMXTR(PRC,XN,ERRFLG) 
IF(ERRFL6,NE,n 60 TO 811 
CALL CLCTR(XN,TRC) 

60 TO 9M 

8fl TYPE 85 • 

85ff FORMAT (24HTR AT INPUTS OUT OF RANGE) 

60 TO Ifll 

9if TR=(TRC*(,8ff MS-,!GMiS<C ))/,8ilf iS 

l«fl RETURN 
END 


NOZZLE PRESSURE RATIO NORMALIZATION 


SUBROUTINE NRMPN 
SCALED FRACTION BRKPT(14) 

DATA BRKPT( I ) ,BRKPT(2) ,BRKPT(3) ,BRKPT(4) ,BRKPT(5) , 

1 BRKPTCG) ,BRKPT(7) ,BRKPT<8) ,BRKPT(9) ,BRKPT(1 i) , 

2 BRKPTd l),BRKPT(12),BRKPT(13),BRKPT(14)/,filf IS, 

3 ,lf||IS,,2MIIS,,3llf IS,,4IIIIS,,5IIIIS,,53685S, 

4 , 611# IS, ,7101 IS,,8I I I IS,, 85 Ills,, 91 Ills,, 95 II IS, 

5 99999S / 

CALL VBNS(BRKPT(I ), BRKPTd 4) ,BRKPT(6>) 

RETURN 

END 


NOZZLE THRUST- VELOCITY CALCULATION, GAMMA=1 .35 


SUBROUTINE CLCFN 
SCALED FRACTION T(14) 

DATA Td),T(2),T(3>,T(4),T(5),T(6),T(7),T<8),T(9), 

1 Tdl),Tdl),Td2),T<13),Td4)/.99999S,.G7l46S,.584l8S, 

2 ,51781s, .45983S, .4I557S, ,385925, .352I9S, .29719S, 

3 .237I9S, .2I312S, .1641 5S,.l 1494S,, III IIS/ 

CALL FNGNKTd)) 


lie 



RETURN 

END 


NOZZLE FLOW CALCULATION AT GAMMA=I,35 


SUBROUTINE CLCWN 
SCALED FRACTION G(1A) 

DATA G(1),G(2),G(3),G(4),G(5),G(6>,G(7),G<8>,G(9), 
! G<1*),G(11),G(12),G(13),G(14)/,676I4S,.G7614S, 

2 .G7614S, ,676 14S, .67614S, ,6761 4S, .6761 4S, .671 I4S, 

3 ,63379S,,55gl8S,.5iiISS»,4216?S,.3f732S,.#IM»S/ 
CALL FNLK1(G(1>) 

RETURN 

END 


ISENTROPIC PRESSURE RATIO NORMALIZATION 


SUBROUTINE NMXTR 

SCALED FRACTION XLOLIM,XHILIM,XINTVL 

DATA XLOLIM,XHILIM,XlNTVL/, 166675,. 9938IS,. 

CALL CBNS(XLOLIM,XHILIM,XlNTVL) 


*16675/ 


RETURN 

END 


I5ENTR0PIC TEMPERATURE RATIO CALCULATION 


SUBROUTINE CLCTR 
SCALED FRACTION T(I4f) 

DATA T d>, T( 2 >, T(3) /, B Bills,. *21 «1 5,. *38725/ 
DATA T(4),T(5) ,T(6) /, 1563! 5, .B7292S, ,*88685/ 
DATA T(7>,T(8) ,T<?) /, 1 1367S, , 1 ! 79PS, . 1 31 695/ 
DATA T(!*).T<1 n.T(12>/.!4484S,.15749S,.16967S/ 


DATA Tn3),T(14>,T(I5)/. 
DATA T(I6),T<17),T(I8)/. 
DATA T(I9>,T<2*),T(2!)/. 
DATA T(22),T(23),T(24)/. 
DATA T(25),T(26),T(27)/. 
DATA T(28),T(29> ,T(3I)/. 
DATA T(3n,T(32),T<33)/, 
DATA T(34),T(35),T(36)/, 
DATA T(37>,T(38),T(39)/, 
DATA T(4*),T(41),T(42)/, 
D^TA T(43),T(44),T(45)/. 
DATA T(46),T(47),T(48>/. 
DATA T(49),T(5I),T(51>/, 
DATA T(52>,T(53),T(54)/, 
DATA T(55),T<56),T(57)/. 
DATA T(58>,T(59),T(6*)/, 


1 81 445.. 192825.. 2*3845/ 

21 4525.. 224895.. 234965/ 

244775. . 254325. .263625/ 

272715.. 28!57S,,29*23S/ 
2987*5,. 316995,, 31 51 15/ 
323*65, ,33*855, .338495/ 

345995.. 353355.. 36*575/ 

367685.. 374665..381525/ 

388275. . 394925. .4* 1 46S/ 
4*79*5,, 41 424S,.42*49S/ 

426655. . 432725 . .438715/ 

444625.. 45*455,. 4562*5/ 

461875.. 467485.. 473* IS/ 

478485.. 483885..489225/ 

494495.. 499715.. 5*4865/ 

519965.. 5151«S,.51999S/ 


117 



DATA T(^:i ) ,T(€S) ,T(63)/.5?>4<52?,.5{>?8fS,,53463S/ 
DATA T(64),T(65),T(66)/,53941S,,544MS,.54883S/ 
DATA T (67), T( 68), T(6<>)/,55347S,.558«7S,, 562625/ 
DATA T(7A),T(71) ,T(72)/,56712S,,57159S,. 576115/ 
DATA T(73),T(74),T(75)/.58f4IS,.58474S,.58?fl5S/ 
DATA T(76) ,T (77 ) ,T( 78)/. 593325, .597755,. 68 1745/ 
DATA T(79),T(8l),T(8I)/.6i59«5,.6| Hi 35 , .61 41 25/ 
DATA T(82) ,T(83) ,T(84) /,61 81 85, .6222115, .6261 <>5/ 
DATA T (85 ),T(86),T (87)/. 6311 55,. 634185,. 657985/ 
DATA T(88),T(89),T(9f)/.641855,.64569S,.6495§5/ 
DATA T(91),T(92),T(93)/. 655285, ,657f35,.66«76S/ 
DATA T (94), T (95 ),T( 96)/. 664465,. 6681 35,. 671 785/ 
DATA T(97),T(98),T(99)/. 675415,. 679015,. 682575/ 
DATA T(1 0i),T(l 01) ,T(1 02)/.686125,.68964S,.693145/ 
DATA T(1 03),T(104),T(1 05)/.696615,.700075,.703505/ 
DATA T(1 06),T(107) ,T( 1 08)/. 706905, .71 0295, .713665/ 
DATA T(109),T(110),T(1 1 1 )/.7 1 7005, .720325, .723625/ 
DATA T(1 12),T(113),T(1 14)/. 726915,. 7301 75,. 73341 5/ 
DATA T (1 1 5), T( 11 6), T( 11 7)/. 736635,. 739835,. 743025/ 
DATA T(1 18),T(119),T(120)/.746185,.749335,.752465/ 
DATA T(121 ), T( 122), T( 123)/. 755575,. 758665,. 761 745/ 
DATA T( 124), T( 125), T( 126)/. 764805,. 767845,. 770875/ 
DATA T(127),T(128) ,T( 1 29)/. 773885, .776875, .779845/ 
DATA T(1 30), T( 131 ),T( 1 32)/. 782805,. 785755,. 788685/ 
DATA T(1 33 ),T(134),T( 135)/. 791595,. 794495,.797375/ 
DATA T(1 36), T( 137), T( 138)/. 800245,. 803095,. 805935/ 
DATA T (1 39), T(140)/. 808765,. 8! 1575/ 

CALL FNGNKT(l)) 

RETURN 

END 


SCALED FRACTION MAPFUN (10,8,12) 


SCALED FRACTION FUNCTION MAPFUN(N,XIN,YIN) 

SCALED FRACTION XIN,YIN,YINCR,XHI,XLO,XFRAC,ZL,ZR,Yl ,Y2, 

1 MAPFUN,XVALS(10,8,12),YVALS(8,12),ZVAlS(lf,8,12) 
C0MM0N/MAPS/XVALS,YVALS,ZVALS,IX(12),JY(12),NX(12),NY(I2) 
NYC- NY ( N) 

NXP=NX(N) 

I=IX(N) 

J=JY(N) 

100 Y1=YIN-YVALS(J,N) 

IF(Y1 .GT.. 000005) 60 TO 110 
IF(Y1 .EQ.. 000005) 60 TO 120 
IF(J.LE.l) 60 TO 500 
JrJ-I 


60 TO 100 

110 Y2rYIN-YVALS(J+l ,N) 

IF(Y2.LT.. 000005) 60 TO 140 
IF(Y2.EQ.. 000005) GO TO 130 
J=J+1 

IF(J,6E.NYC) GO TO 500 
Y1 = Y2 


GO TO 110 
120 YINCRr. 000005 
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GO TO 2M 
!3i YINCRr.999995 
GO TO 2if 

14« YIHCR=YI/(Y1-Y2) 

2ff XLO=XVALS(I,J,M)+YINCR*(XVALS(I,J+l,N)-XVALS(I,J,M>) 
IF(XIN.GT.XLO) GO TO 22# 


22 # 


25# 

3 ## 


4 ## 


9#l 


IF(XIN.EQ.XLO) GO TO 23# 

IFd.LE.l) GO TO 5## 

IrI-1 
GO TO 2## 

XHI=XVALS(I+1 ,J,N)+YINCR*(XVALS<I+I , J+I ,N)-XVALS(I-H , J,N>) 
IF(XIN.LT.XHI) GO TO 3## 


1 = 1+1 

IF(XIM.EQ.XHI) GO TO 23# 

IFd.GE.NXP) SO TO 5## 

XLOrXHI 
GO TO 22# 

MAPFUN=ZVALSd,J,N)+YINCR*<ZVALSd,J+l,N)-ZVALS(I,J,N)> 

GO TO 4## 

ZR=ZVALSd+l ,J,N)+YINCR+(ZVALS<I+I ,J+1 ,N)-ZVALS(I+I , J,N)> 

ZL=ZVALSd,J,N)+YINCR*(ZVALSd,J+l,N)-ZVALSd,J,N)) 

XFRAC=(XIN-XL0)/(XHI-XL0> 

WAPFUN=ZL+XFRAC*(ZR-ZL) 

IX(N>=I 
JY(m=J 
R ETUR N 

CALL MOOR(N,XIN,YIN) 

RETURN 

END 


MAP OUT OF RANGE ROUTINE (MOOR) 


SUBROUTINE MOOR (N, XI N, YIN) 

SCALED FRACTION XIN,YIN 

C TEST FOR OPERATE MODE 

CALL QRAMidLOC) 

IF(IL0C,EQ.4) CALL QSHdERR) 

TYPE «##,N,XIN,YIN 

Si# F0RMAT(/7HMAP NO.,I3,2#H INPUTS OUT OF RANGE/SHXIN = ,S7, 
I 8H YIN s ,S7/) 

PAUSE 

IFdL0C,EQ,4) CALL QSOP(IERR) 

RETURN 

END 
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APPENDIX F 


SUMMARY OF SCALED EQUATIONS SOLVED ON ANALOG COMPUTER 


W21 = C(1)*DW21 dt' + C(2) 


(FI) 


T21 = r C(3)* ^ 

JO W21 


(F2) 


P21 = 10 *C(26)*W21*T21 


(F3) 


W3 = ^ 10 *C(11)*DW3 dt’ + C(9) 


(F4) 


rV 

Tg = 10 *C(12)* dt’ + C(10) 

JO W3 


(F5) 


P3 = 10 *C(34)*W3*T3 


(F6) 


fV 

W4 = 10 *C(15)*DW4 dt’ + C(13) 


(F7) 


= f ^ 10 *C(16)* df + C(14) 

Jn W4 


(F8) 



P4 = 10 *C(38)*W4*T4 


ft’ 

W5 = / 10*C(19)*DW5 dt’ + C(17) 

Jo 


T5= 10 *C(20)* dt' + C(18) 

Jo W5 


P5 = 10 *C(41)*W5*T5 


W6 = p' 10 *C(23)*DW6 dt' + C(21) 
^0 


T„ = f 10 *C(24)* dt’ + C(22) 

® Jn W6 


P6 = 10 *C(43)*W6*T6 


W7 = r C(33)*DW7 dt' + C(35) 
Jo 


T7 = r C(48)* dt' + C(49) 

Jo W7 


P7 = 10 *C(36)*W7*T7 


(F9) 


(FIO) 


(Fll) 


(FI 2) 


(F13) 


(F14) 


(F15) 


(F16) 


(F17) 


(F18) 
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W32 



C(29)*DW32 dt'+C(3l) ) 


ft' 

132 = / C(46)* T32m ^ 

JO W32 

P32 = 10 *C(32)*W32*T32 


W35 



C(25)*DW35 dt'+C(27) 


T35 


= f C(44)* T35m +c(45) 
JO W35 

P35 = 10 *C(28)*W35*T35 


W22 



DW22 dt’+C(6) 


T22=f^ C(7)*T2^ dt-+C(8) 
JO W22 

P22 = 10 *C(30)*W22*T22 


PCNTl 



C(37)*DTQT1 dt'+C(39) 


(F19) 

(F20) 

(F21) 

(F22) 

(F23) 

(F24) 

(F25) 

(F26) 

fF27) 

(F28) 
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ft' 

PCNT2 = / C(40)*DTQT2 dt' + C(42) 

Jo 

(F29) 

ANl = C(56) 

(F30) 

AN2 = C(54) 

(F31) 

P2 = C(55) 

(F32) 

WDF = C(57) + Output of controller 

(F33) 


If SUPER="TRUE", 

P21 = 10*C(52)* [C(51) + 10*C(50)*C(53)*PCNT2] (F34) 
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APPENDIX G 


COMPUTER LISTING OF TURBINE DATA CONVERSION PROGRAM 


C ION CHT( ]J ) T ( c ) , k»R 7P ( 1G» 9) ,h TA( 10 »9 ) , Fh ( IC,9 ) 

C PP l( 12) ,V*RTP2( 12) ,F1 (12) ,F2(12 ) 

C )FFNS I1N RATinC 1C,9) ,nHTi (1C,9) ,V<RTP1 (10,9) ,Ll,RT (10,9 ), 

I HPM10,9),NR1F1( 10,9),NR7F2(IC,9),LI"PTP2(10,9) 

REAL NRT,l kR7,Nfi 7Fl,NRTF2,LViRTF2 

^ AREl 1ST / FEA7/DH7 / SPFEC/F:R7 /FI CV\/V»R7P/EF F/£T4 

NAVEL I ST /PAR AM 1/ PR 1 / PA R AM2 / )>R 7P2 

CP = . 28595 

GAV= 1 .33 12 

R£Ar(5,HEA7 ) 

PEAC(5,SPEEn ) 

Rf AP(5,PLCM 
PEAC(5,cFF ) 

R 1 AC( 5,PARAM 1 ) 

R f AC ( 5,PAR AM 2 ) 

C 

c PRINT TIJReiNE DATA TABLE HEAOINC: 

C 

(nP ITF( t, ICC ) 

IOC FCPMAT( IFI ,50X,18HTURBINE DATA TABLE) 

WP 17E ( 6, 1C 1 ) 

101 FCFMAT(1FK,1AX,3HDFT//7X, 3HNR T,5X,3HETA/15X,4FtV*RTP/15X,2HPP) 
)»R 1TE(6, 1C2) (DHTd) ,1=1,10 

102 FCPMAT( IFK ,14X,1C(F6.A,3X) ) 

C 

c CALCLLATF TURBINE PRESSURE RATIC 

C 

CC 2 J=1,S 

tC 1 1=1,10 

C)CF=1./1.81SC£ 

ETTCF=.972AA 

PP(I,J)=(1.~DHT(I)/(CP*ETA(I,J)*DHCF*ETTCF))**(GAI>'/(GAF-1.)) 

1 CCKTINUF 
C 

c PRINT TUReiNE DATA TABLE 

C 

WRITE(6, ICE) NRT(J),(ETA(I,J),l = i,iO) 

102 FCRMAK 1HK,6X,F5.1,3X,1C(F6.4,3X)) 

WPITE(6,lCO URTPd ,J) ,I =1 ,10) , ( PR ( I , J ) , 1=1 , 10 ) 

10 A FCRMATI IH J , 1 4X , I C ( F 6 . 3 ,3 X ) / 1 5 X , I 0 ( F6 . 4 ,3 X ) ) 

2 CCMINUF 
C 

C PRINT TUREINE PRESSURE RATIC MAP HEADING 


C 

C 


WP ]TE( 6, 1C 5 ) 

105 FCPMATI IFl ,45X ,26HTUR6 INt PRESSURE RATIO MAP) 
WP 1TE(6, ICC ) 


IOC FCPMAT(IFK,14X,2HPR//7X,3HNRT,5X , 3HDH T/1 5 X ,4HW FT P/ 15X ,4HLWRT/ 15X, 
1 <FWTPN/1AX, 5HNR TFl/13X,6HLWRTF2/14X,5HNRTF2) 

WRlTE(C,in) (PRKK) ,K = 1,12) 


INTERPOLATE TABLE FCR CONSTANT PRESSURE RATIOS 
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r 


cr f j=i,<; 

cr ^ K= 1 , 

PP 2 = PR L(K ) 

F ](K )= SQPK l.C-PP 2»J*PP 2) 

F2(»^)= SOJ^T/l.C- f"Rl(K) ((G4W -1.0) / G4^)) 

I F ( ( PR l( K ) .r,T, ^K( I, J) ).rR. (PKKK) .LT.PR (10, J) )) CC TC 5 
rr 3 1= 1, 1C 

IF (PR KK HPP ( I , J > ) 3,4,4 
2 CCN T IMUE 

4 RMIC(K,J)=: (PF ( I- 1, J )~PR 1 ( K) ) /( PR (I -1 ,J ) - FP ( I , J ) ) 

Cn 1 (K, J )= DHT( I- n^RATIO (K , J) « (DH T( I) -r^HT (I-l ) ) 

WR^Pl(K,J)r WRTP(I-l,J) + PATIf(K,J) + ( R T ( I , I ) - W RT P ( I - 1 , J ) ) 
LVP7(K,J)= 0HT1(K,J) / ^RT(J) 

U1PM(K,j)r WRTPKK,)) / KRl(J) 

NF1FI(K,J)= MRT(J) / F1(K) 

NR1F2(K,J)= NP7(J) / »^2(K) 

LVRTF2(K,J)= LWRT(K,J) / F2(K) 

5 CCN^IMJF 

WF ITF(6, 10?) NRT(J), (0HT1(K,J) ,K = 1 ,12) 

107 FCPVAT(lhK,6>,F5.1,3X,12(F6.4,3V)) 

UPTTF(6,lll) (WRTPI(K,J), K^l^lZ) , (L^R^(K,J), »<=1,12), 

1 (WTPN(K,J), K = 1,12), (NR1FI(K,J), K = l,12), ( LWPT F2 ( K, J ) , 

2 K = 1,12), (NPfF2(K,J), K=1,12) 

111 format (IFJ, 14X, 12(F^.2 ,3X) / I'SX, 12(t^.3) / 15X , 12(F6.4, 

1 :>) f 15X, 12fF^.l, 3>) / 15X, 12(^7.6, 2>) / L5X, 12(F6,1, 3X)) 

f CCNTINUF 
C 

C PRINT TUPeiNE FLGVs MAP HEADING 

C 

V^R ITE(6, ICF) 

ice FCRMAT( IM ,45X,16HTU«8INL FLOU MAP) 

WR 1TE( 6, IGF ) 

lOS FCRMAK IFK , 14X , 4HWRTP / /7X , 3HNRT t 5X ,3HDHT ) 

WFITE(6,11C) (VsR.TP2(X), K = l,I2) 
no FCPMAT( 1HK,14X,12(F6.3,3X) > 

C 

C INTERPOLATE TABLE FCR CONSTANT FLC^ 

C 

cr 12 J=l,9 

CC 12 K=1,12 

CFT1(K,J)= C.C 
12 CENT INUE 

CC 11 J^ltS 
CC 1C K=l,12 

I F ( ( ^RTP 2( K) .LT, WR TP (1 , J) ) . CR. ( WRTP2 (X)* GT .WPT F( 12 , J) ) ) GC TC 10 
cr 7 1=1, 1C 

IF(VsRTP( I, J )-VsRTP2(K) ) 7,8, E 

7 CONTINUE 

8 1F(VPTPC I- I, J)-\^RTP( I , J) ) 9,10,0 

9 P AT 10(K, J )= { WRTP ( I-' 1 , J)-t^RTP2 (K) ) / ( WRJP ( I - I , J )--W P T P ( I , i) ) 

CM 1(K, J )= DHT (I - DORATI C ( K,J) >«^(DHT( I l-DHT C I- 1 ) ) 

IC CCNT INUE 

VsR 1TE( E, 1C7 ) NRT(J), ( CH T 1 ( K , J ) , K = i , 12 ) 

11 CCNT INUE 
STCP 
E^ C 


125 



rMU)= •Clt.C?, .C3fC35,.C4t.C483,.0bS 
^SPF^C NkK 1)=4 1. 3 6,49.n,f 7.9Ct66. 17,74,4^, 
l^,qP,P2.n,<;C.<;q,cq.26$ 

$FLCW ^RTP( 1, l) = 65.,8P.,lC3^,lC7,,llu*2, 

1 13 1 1 4. 6 , I I 5. 5 tl 16* ,1 16.3 , 

^P1P(1,2) =6 4.0,64. ,<;8.,1C3.,LC 7., 

I 1C .6, 1 12. , 113. ,114. , 11 5. , 

^FTF( 1,3) =7 ".,6 2.,^5. ,1CC. ,1C3.5, 

1C8., lie. 1,111.5,112.8,113.3, 

WP1P(1,4) =6^.,82.,52. ,^8. ,1C1. , 

1C5.6,1C8.,1CS.2,11U7,112.5, 

WP1P(1,5) =65. ,82. ,91. ,65. ,5€. 5, 

ia3.,lC6.5,lC7.5,llC.,111.5, 

WF1PC1,6) =7 1 82. ,9C. 5 ,54. ,97. , 

1C2.23,1O5.,1C7.,1C0.,111., 

^P1P(1,7) =71. ^,82. ,50.5,53.5,56.5, 

IC2., IC4. , 1C6. , 1C9. ,IIC.5 , 

KRTP(1,8) =75.,82.5,5C.,93.,96., 

ICC., 1C2.5,1C4.,107.,IG9. , 

UF1P(1,Q) =76. ,83. ,5 C. ,52. 5 ,95. , 

95., 102. , IC3. ,1C6. ,107.51 
iEFF El/S(l,l)= .88, .86 ,. 76 , .71 ,.64, 

. 55 ,. 55 ,.5^ ,. 55 ,. 55, 

El/Ml, 2 ) = .82, . 89 ,. 87, .84, .81, 

. 75 , . 7 1 ,. 67 ,. 61 ,. 55 , 

El^(l,3) = .78,. 8 8 3, . 5C3 ,. 5 , . 887 , 

. 6 55 , . 63 , . 6 1 5 ,. 78 ,. 75 , 

E7^(l,4) = .64,.84,.5C,.51,.913, 

.9C3,. 852 ,. 8 E2 ,. 855 , . 835 , 

El^(l, 5 ) = . 55 , . 80,. 877 ,. 855,. 907 , 

. 922, .9 1 6 ,. 5 1 1 , . 501 ,.P9, 

E17(l, 6 ) = . 52 , . 7 5,. 645 , . 88 , . 897 , 

. 518, . 525 , . 522 ,. 513 ,. 905, 
ei/i(l, 7 ) = .50 , . 73 ,. 64 ,. 865 , . 89, 

. 5 1 , . 925 , . 525 , . 5 13 , . 91 , 

ET/i(l,8) = . 45, . 7, . 8 ,. E35 , . 86 , 

. 853, .9C7 ,. 522,. 53 ,. 52 , 

El^(l,9) = . 45 ,. 65 ,. 76 , . 8, . 825 , 

.86 5, . 8 52 ,. 5C2, . 52 5 ,. 92 5i 

$PAPA^'J PP1(1) = .1,.2,.3,.4,.5,.6,.7,. 75 ,.8,.85,.9,.95^ 
fiPAPAM; ^RTP2(1) = 50.,6C.,7C.,8C.,90.,95., 

10C.,102.2,lC5.,lC7.5,llG.,115.t 
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APPENDIX H 


TURBOFAN FUEL CONTROL SIMULATION 


In order to investigate the dynamics of turbojet and turbofan engines, it is often 
necessary to simulate the engine fuel control system, as well as the engine. The use of 
the HYDES program dictates that the fuel control be simulated on the analog computer. 

A portion of a proposed fuel control was simulated on the analog computer to allow 
transient variations in the operating point of the selected turbofan engine. The primary 


function of the control is to maintain the low- pressure- turbine rotor speed N^2 at a 


demanded value N^g^ ^lem' 
of the speed control loop. 


Figure 16 illustrates, in block diagram form, the operation 


The measured speed N^2 m subtracted from the demanded value, which is 
scheduled as a function of power lever angle a and inlet temperature T2 The error 
is amplified and then reduced by the output of a derivative-plus-lag circuit whose input 



I From acceleration 
schedule error 


Figure 16. - Turbofan low-pressure-tgrbine rotor-speed control. 



Figure 17. - Turbofan acceleration control. 
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is fuel valve position y^. The feedback provides proportional-plus-integral dynamics 
between the valve position and the speed error. For near-steady-state conditions, the 
MIN circuit (output equals smallest of its inputs) transmits the speed error signal to 
the fuel valve servomechanism. 

In order to combine fast thrust increases with overtemperature and compressor 
stall protection, an acceleration schedule feature is included in the control and is illus- 
trated in figure 17. 

A sudden increase in the power lever angle a causes a corresponding increase in 
the speed error signal. The MIN circuit then transmits the acceleration error signal to 
the fuel valve servomechanism. An acceleration schedule parameter (p is scheduled as 
a function of a , high-pressure-turbine rotor speed N^.j and compressor inlet 
temperature T^ . . The output of the acceleration schedule is multiplied by compres- 

sor discharge pressure to determine the acceleration fuel flow and corresponding valve 
position. The amplified position error is then transmitted to the fuel valve servomech- 
anism, giving a first-order lag response of valve position to the command acceleration 
position. The engine then accelerates along the acceleration schedule until the speed 
error signal becomes smaller than the acceleration position error signal. 

Riysical limits are built into the valve position. That is y^. ^yp <yp j^ax’ 

The speed control feedback signal (fig. 16) is also limited to y^g ^ ypg <ypB max' 

For the turbofan transient demonstration, the low-pressure-turbine rotor speed 
demand N |,2 was ramped from its initial value to the final value in 50 milliseconds 
(5 sec computer time). The acceleration schedule was fit by a function of compressor 
inlet temperature and speed. The following equations were implemented on the analog 
computer: 




'^2. 1, m ■ *^2.1 ■ 


^ %,i 


(HI) 




(H2) 


*^2. 1, m^ ^ "^2. 1, i 


(H3) 
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!.m=50| (Pj - Ps,^)dt . ^ 


(H4) 


Wp = 4.653 (yp + 0.0846)' 


(H5) 


= 1. 664 P 

Jo 


MIN 


^t2 


accel 


dt + yp j 


(H6) 


t2 


25.91 


N, 


t2, des 


(^t2,dem “ ^t2, m) “ ^FB 


(H7) 


FB 


= 18. 014 r 
Jo 


Yf - -^)dt 
, ^ 9. 007 / 


(H8) 


'accel F, accel “ ^F^ 


(H9) 


^F, accel 


w 


F, accel 


4.653 


1/2 


(HIO) 


.r, 

W-r-, = 

F, accel 


3600 


(Hll) 
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<p 


= K 


N. 


<P 


+ 33 


tl, m 


N. 


tl, des 


^std 

2. 1, m 


/o. 3124 + 0.6895 

^std 


(H12) 


K 


<p , nom 


-5.4 
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